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Abstract. This work is devoted to the study of a family of almost periodic 
one-dimensional Schrodinger equations. We define a monodromy matrix for this 
family. We study the asymptotic behavior of this matrix in the adiabatic case. 
Therefore, we develop a complex WKB method for adiabatic perturbations of 
periodic Schrodinger equations. At last, the study of the monodromy matrix 
enables us to get some spectral results for the initial family of almost periodic 
equations. 

Resume. Ce travail est consacre a I'etude d'une famille d'equations de Schrodin- 
ger quasi-periodiques en dimension 1. Nous definissons une matrice de mo- 
nodromie pour cette famille. Nous etudions le comportcmcnt asymptotiquc do 
cette matrice dans le cas adiabatique. A cette fin, nous developpons une ver- 
sion de la methode WKB complexe pour I'equation de Schrodinger periodique 
avec une perturbation adiabatique. Enfin, I'etude de la matrice de monodromie 
nous permet d'obtenir des resultats spectraux pour la famille d'equations quasi- 
periodiques initiale. 



1. Introduction and statement of the main results 
In this paper, we study the family of Schrodinger equations 

(1.1) — — ^(x) + (V (x -(!)) + W (ex)) i)(x) = E^( x), x e M. 

dx'^ 

Here 

• e M is a parameter numbering the equations, 

• V and W are two periodic functions satisfying 

(1.2) V{x + l) = V{x), W{i + 2Ti) = W{0, x,ie% 

• £ is a fixed positive number. 

The ratio of the periods of the potential in (1.1) is equal to 27: /e. If 27r/£ Q, 
the potential V [x — (j)) -\-W {ex) is almost periodic. 

Our ultimate goal is to understand the spectral properties of the family (1.1) 
assuming that V is a real- valued measurable bounded function of a: e M and that 
W is analytic in ^ e C. 
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It is well known that for one-dimensional periodic Schrodinger equation, the key 
to the understanding of its spectral properties is the monodromy matrix. It is not 
possible to define this object for a single almost periodic differential equation. In 
this paper, we will show that, for a family of almost periodic differential equa- 
tions of type (1.1), a monodromy matrix can be defined and used for its spectral 
analysis. 

Of course, a thorough understanding of the structure of monodromy matrix 
is crucial. Therefore, in this paper, we will need to take the adiabatic limit e ^ 0. 
As we shall see, the monodromy matrix always enables us to reduce the study of 
(1.1) to the study of a difference equation. In the adiabatic case, this equation 
takes a standard form and can be effectively analyzed. 

1.1. Reduction to a difference equation. We begin by introducing several 
natural notions related to the family of equations (1.1). 

1.1.1. For any fixed 0, let ijji^2{x, ^) be two linearly independent solutions of 
equation (1.1). We say that they form a consistent basis if 

(1.3) ^^Jl^2{x, + = i>i,2{x, (I)), Va;,(/», 
and if their Wronskian defined by 

(1.4) (^1, ^2) = -0^-02 - '02'01- 

is independent of (/). For convenience, we shall assume that 

(1.5) w(:^,, ii'2) = i, v</.. 

Remark that the solutions verifying the Cauchy conditions 

^i(0,</.) = 0, ^i(o,</.) = i, 

form a consistent basis. However, as we shall see in the next sections, this choice 
of a consistent basis is not the most natural, and, for a while, we do not choose 
any particular construction of the solutions ipi and ■02. 

1.1.2. In view of (1.2), the functions tpi^2{x+27r/£, 0+27r/£) are also two solutions 
of equation (1.1). Therefore, one can write 

(1.6) * (x + 27r/e, + 27r/e) = M (0) * (x, 0), 
where ^ is the vector 

(1.7) ^^{x,(|>) = {^P,{x,(|>),^P2{x,(^>)), 

where ^ is the symbol of transposition, and M (0) is a matrix with coefficients 
independent of x. We call this matrix the monodromy matrix corresponding to 
the consistent basis tpi^2- 

Note that this is a natural generalization of the canonical definition of the mon- 
odromy matrix: if, in the Schrodinger equation (1.1), V = 0, then the potential 
is periodic, the basis solutions -01^2 are independent of 0, and formula (1.6) de- 
fines the standard monodromy matrix corresponding to two linearly independent 
solutions of the periodic Schrodinger equation, see [21]. 
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It follows from the definition that, for any consistent basis, 

det M ((/)) = 1, M((/)+ 1) = M (</)), V(/). 

1.1.3. Set 

(1.8) /i=— modi. 

s 

Let M be a monodromy matrix corresponding to a consistent basis ■0i^2- In Sec- 
tion 6, we shall see that the spectral analysis of (1.1) can be reduced to the 
investigation of the solutions of the equation 

(1.9) x{'P + h) = M{cl))x{4>). VxeM. 
Assume that x is a fundamental solution i.e. that 

(1.10) detx(</') = l, V(/)eE. 
One can easily obtain that 

(1.11) * (x + 2-Kn/e, (l)) = x X'^{(l> - nh) * (x, ^/i) - nh). 

It is well known that the spectral properties of one dimensional Schrodinger equa- 
tions can be completely described in terms of the behavior of its solutions as 
X ±oo (see [12]). Remembering that is 1-periodic in (f), we see that rela- 
tion (1.11) reduces the analysis of the behavior for x — )■ ±oo of the solutions of 
the input differential equation (1.1) to the analysis of solutions of the difference 
equation (1.9). 

The passage from equation (1.1) to (1.9) is done in the same spirit as the mon- 
odromization transformation [8]. For example, this latter transformation reduces 
the spectral analysis of a difference Schrodinger equation to the analysis of a differ- 
ence equation of the form (1.9); in this case too, the matrix defining the resulting 
equation is the monodromy matrix corresponding to one of the fundamental so- 
lutions of the initial difference equation. 

Here, as in the case of the monodromization, one can define Bloch solutions for 
the almost periodic family (1.1), and show that these Bloch solutions can be 
constructed in terms of Bloch solutions of the difference equation (1.9) (for the 
definition of the latter, see [8]). These ideas proved to be eff'ective for the spectral 
analysis of the Harper equation [6]. We will discuss and use these ideas in section 6. 

1.2. Complex WKB method. To make the above ideas work, one has to find an 
effective way to control the monodromy matrices for family of equations (1.1). For 
this, it suffices to have a good enough control of some consistent basis solutions. 
In this paper, we construct and study certain consistent basis for small values 
of £ using the ideas of the complex WKB method. The complex constructions 
appear to be inevitable since some of the coefficients of the monodromy matrices 
are exponentially small when s ^ 0. The results and the constructions developed 
first for (1.1) are very general, and hopefully, they may be applied to study a wide 
class of problems for adiabatically perturbed one-dimensional periodic Schrodinger 
equations. So, we formulate the results in a more or less general form without 
assuming that the function W is periodic. We also note that our assumptions on 
the regularity of V and W can be essentially weakened. 
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1.2.1. We begin with some preliminary remarks concerning the standard version 
of the complex WKB method developed for ordinary differential equations of the 
form (see, for example [11]) 

(1.12) -e''—ij{z)^U{z)4!{z) = Eij{z), zeC. 

Here, U is an analytic function of z, and e is the standard small semi-classical 
parameter. For (1.12), the complex WKB method consists in the description 
of certain canonical domains in the complex plane 2 G C, and in constructing 
analytic solutions of (1.12) having standard asymptotic behavior for e — > on 
a given canonical domain. For a given canonical domain K, one constructs two 
linearly independent solutions il)± with the asymptotics 

Here Zq is a fixed point in K and p is the complex momentum defined by the 
relation 

(z) + U{z) = E 

(i.e. by the symbol of the equation (1.12)). The complex momentum is the 
main analytic object of the WKB constructions. The canonical domain can be 
characterized as a domain such that, through any point of it, one can draw a 
smooth curve along which the function Im J^^ p {z) dz is monotonically increasing. 
For precise definitions, we again refer to [11]. 

1.2.2. For the family of equations (1.1) the role of the complex momentum is 
played by the function k [^p) satisfying the relation 

(1.13) 8{k) + W {^) = E, 

where £ (k) is the dispersion law of the periodic Schrodinger equation 

(1.14) -'ij;(x) + V (x) = Eiplx), xeR. 

du^ 

In other words, 

(1.15) K{^) = k{E-W{^)), 

where k (E) is the Bloch quasi-momentum of (1-14), see subsection 2.2. As p{z), 
the complex momentum k is a multi- valued analytic function. Its branching points 
are related to the branching points of the quasi-momentum. 

1.2.3. The simplest asymptotic behavior of analytic solutions of (1.1) is predicted 
by the already classical ansatz suggested by V. S. Buslaev. In his papers (see, 
for example, [1], [2], [3], [4]), V. S. Buslaev has studied the asymptotic behavior 
along the real line x G M of solutions of the equation 

(1.16) — —il)(x) + U(x,ex)il)(x) = Eil)(x), xeR. 

dx^ 

In this equation, U (x,^) is a sufficiently smooth function periodic in x. Between 
the branching points of the complex momentum, the leading term of the asymp- 
totics of two linearly independent solutions has the form 

±- r"" K (^) d^ 

(1.17) il)± A±{sx) e E Jexo ^''^ ^p±{x,ex), e ^ 0. 



Here, p±{x,^) are the periodic component of the Bloch solutions (see subsection 
3.1 and 3.3) of the equation 

(1.18) -^V' (x) + U {x, (x) = Ei;{x), xe R, 

where the parameter ^ G M is frozen and ^±(0 two "normahzation constants" 
described in terms of the periodic components and the complex momentum. 

Buslaev's ansatz played a very stimulating role for our constructions; moreover 
the asymptotic formulas of the complex WKB method developed here look similar 
to the ones of Buslaev. However, Buslaev's constructions can not be used as a 
base for the constructions of the complex WKB method needed here. 

1.2.4. We will now formulate the main theorem of the complex WKB method 
for the family of equations (1.1). For this rewrite (1.1) in terms of the variables 

U = X — (p, ip = £(j). 

The equation takes the form: 

(1.19) -—'ip{u) + {V{u) + W{£u + if))^{u) = E^{u), liGR, 

and the consistency condition (1.3) comes to 

(1.20) iii^2(u - I, LP + £) = ibi^iiu^if), ip,u&'R. 

Recall that the function W is analytic. So, we may choose the parameter -p to be 
complex. The goal of our version of the complex WKB method is twofold: 

1. describe certain canonical domains on the complex plane e C, 

2. on any given canonical domain, construct a consistent basis of solutions of 
(1.19) having the simplest asymptotic behavior when e — > (li being in a 
fixed interval [-[/, [/] C M). 

When describing the constructions of the complex WKB method we shall refer to 
the equation 

(f 

(1.21) -—il) + V{u)il) = 8il), mgM, 
in which 

(1.22) £ = E-W{^), 

and (y9 e C is fixed. The main theorem of the complex WKB method in the 
adiabatic case is 

Theorem 1.1. Fix [/ e M. Let K be a bounded admissible canonical domain 
for the family of equations (l-l). For sufficiently small £, there exist a consistent 
basis /-t defined for « e M and (p G K and having the following properties: 

• for any fixed u eM., the functions f±{u, ip) are analytic in (p E K. 

• for —U <u<U, the functions f±{u,(p) have the asymptotic representations 

(1.23) f±{u,Lp) = e'^-s Io'^'^^{^±{u,Lp) + o{l)), e^O. 

Here ^-t are the canonical Bloch solutions of equation (1-21) corresponding 
to the domain K. 
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• the error estimates in (1.23) are uniform in and u and may differentiated 
once in u. 



Roughly, a canonical domain is a domain such that through any point of it, one 
can draw a smooth curve along which the function Im JJ'^ k dip is monotonically 
increasing and the function Im J^^ {k — tt) d(p is monotonically decreasing. An 
admissible sub-domain of a canonical domain is this canonical domain without a 
fixed (5-vicinity of its boundary. The Bloch solutions mentioned in the theorem 
have the form 

Here ip± are two linearly independent Bloch solutions of equation (1.21) defined 
for if ^ K. The factors A± are constructed in terms of the periodic components 
of these Bloch solutions by the same formulae as in the case of (1.17). The reader 
will find precise definitions in section 4. 

We prove Theorem 1.1 in section 3. Let us emphasize here the crucial role played 
by the consistency condition; indeed, it is this condition that allows us to control 
the dependence of the solutions on ip and leads to the asymptotics (1.23) and to 
the definition of the canonical domains. 

In section 4, we consider the equation 

d^ 

(1.24) — — il)ix) + V ix - (I)) iex)il)ix) = Ell) ix), a; G M, 

dx'^ 

analyze the analytic properties of the complex momentum n and describe all 
the related geometric constructions and some of the canonical domains. These 
canonical domains are shown on Fig. 7 and Fig. 9. They remind us of the 
canonical domains arising in the semi-classical analysis of Harper equation 

/-, r^rX i'i^ + h) + tp(x - h) , , , , , / N 

(1.25) — ' ^ + \co&xiIj{x) = EiIj{x), 

in which h is the semi-classical parameter (the case of A = 1 was investigated in 
[7]). We will discuss this similarity later. 

1.3. Asymptotics of the monodromy matrices. For any fixed m G M, the 
WKB solutions described in the Theorem 1.1 appear to be analytic in a whole 
strip 

of the complex plane. Indeed, being analytic in a given canonical domain, they 
can continued analytically outside just by means of formulae (1.20). As a result, 
one can see that the corresponding monodromy matrices are analytic in the same 
strip. In sections 5 - 6, we apply the above constructions to equation (1.24). The 
structure of the canonical domains for this equation and, thus, the asymptotics of 
the monodromy matrix strongly depends on the value of the spectral parameter 
E. Denote by E\ and E^ the ends of the first spectral interval of the periodic 
Schrodinger equation (1.14). For the sake of the definiteness, we assume that 

(1.26) E-l<Ei-5, El + 5 < E + 1< ^2 - 5, E2<Ei 

where ^ is a fixed positive number. Starting with a consistent basis corresponding 
to one of the canonical domains, we prove in section 6 the 
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Theorem 1.2. Under the condition (1.26), for sufficiently small positive e, there 
exists a consistent basis and a positive number Y such that the corresponding 
monodromy matrix is analytic in the strip \lm(p\ < Y and can be represented in 
the form 

do ^0 + ^le^^ 

(1.27) M = T - \ I -T-K 



Co + cie do + die^^ + d_ie 



Here 



• T is a constant 2x2 diagonal matrix. 

• z = 27r(j) + C (C is a real constant), 

• the coefficients a, b, c and d admit the asymptotics: 

do = f C0S((/)1 +o(l)) (l+o(l)) , = -^(l+o(l)), d_i = -^(1+0(1)), 

bo = ie^^^ (1+0(1)) , bi = -iti (i+o(i)) , 
Co = -ie^^^ (1+0 (1)) , ci = iti (1+0 (1)) , 
ao = te^^i (1+0(1)) , 

where 0i, t and ti are nonzero real coefficients depending only on E. The asymp- 
totics are uniform in E and in ip. 

Let us discuss the above result. First, we describe the coefficients cj)i, t and 
ti. The discussion given below easily follows from formulae (5.1). Remember that 
the complex momentum k ((p) is a multi-valued function. Under the condition 
(1.26), it has branching points at the real line (see Fig. 2-4). They are described 
by the relation 

El + cos = E. 

There is only one branching point on the interval (0,7r). Denote it by cpi. All the 
other real branching points are the points (p = ±(pi + 27r/, I G Z (see Fig. 2-4). 
There is a branch of the complex momentum that belongs to iR^ on the intervals 
(— + 27rl, (pi + 27r/) and to R_|_ on the intervals {ipi + 27rl, 27r — (^i + 27r/), / G Z. 
Denote this branch by k*. 

Consider the real isoenergy curve (1.13) (i.e. the intersection of the curve (1.13) 
and the plane {p,k) E M'^). Under the condition (1.26), it corresponds to the 
figure given bellow. This curve is 27r-periodic in p and n. The branch of this 
curve connecting the points pi and 27r — ipi is described by the formula 

K = K^{(p), Pi < P < 277 - pi. 

The coefficient t is the tunneling coefficient describing the tunneling between the 
branches of the isoenergy curve separated by the points —pi and pi. The co- 
efficient ti is the tunneling coefficient describing the tunneling of between the 
branches of the isoenergy curve separated by the points Ki and 27r — Ki. The 0i 
is the phase integral associated to the closed branch connecting the points p)i and 
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2n — (pi. More precisely, 



(1.28) 

(j)i = - / K,*dip, t = exp{- / K^dip), = exp(- / {ip* — 7r)dK). 

Here. (P*{k) is a complex branch of the isoenergy curve connecting the points 
(fi = 71, K = Kl and (/3 = TT, K = 2% — Kl. Along this branch (P*{k) — tt G iM^. Note 
that the tunneling coefficients t and ti are exponentially small in and that 
the phase integral 0i is real. 

The integrals defining ^i, ti and t can also be written as contour integrals in 
{k, (f) -complex plane. For (pi (resp. t, ti), one integrates the form Kdcp along the 
contour 71 (resp. 72, 73). The contour 72 (resp. 73) is imaginary in k (resp. in 
ip — tt) and real in ip (resp. in k) (see figure below). 




The topological structure of the real isoenergy curve (1.13) is the same as the 
one of the curve cosp + cosx = E with —2<E<2. The latter one is the 
real isoenergy curve defined by the symbol of Harper equation (1.25) (A = 1). 
This leads to a similarity of the above mentioned canonical domains for equations 
(1.25) and (1.1) in a neighborhood of the real line. 

We prove Theorem 1.2 in section 5. The calculations leading to Theorem 1.2 use 
standard ideas of the classical complex WKB method. Due to the similarity of 
the canonical domains, they are parallel to the calculations of made in [5] for the 
semi-classical analysis of monodromy matrices for Harper equation. 

Let us note that it is Theorem 1.1 which reduces the asymptotic analysis of the 
almost periodic family of equations (1.1) to the investigation of objects, most of 
which well known in classical versions of the complex WKB method. 

We finish this subsection by discussing the asymptotic analytic structure of the 
monodromy matrix (1.27). In the Theorem 1.2, we have described the asymptotics 
of the first Fourier coefficients of the monodromy matrix. We can show that if the 
periodic Schrodinger operator with the potential V has infinite number of gaps 
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in its spectrum, then the Fourier series are also infinite. However, in the strip 
|Im (^1 < Y, the leading term of this matrix is a trigonometric polynomial of order 



one. 



Consider a family of monodromy matrices generated by Harper equation in course 
of the monodromization renormalization procedure, see [9]. The matrices of this 
family have the form 



(1.29) 



2A cos 2; s + te 



te^^ st/X 



where a, s, t are constant parameters. Since the determinant of this matrix equals 
to one, only two of them are independent. Comparing (1.27) with (1.29), we see 
that one can hope to carry out the spectral analysis of the family of equations 
(1.1) by means of the monodromization method developed for Harper equation. 

1.4. First spectral results and main conjectures. 

1.4.1. Let Js be the interval of energies E defined by (1.26). In Section 6 we 
show that Theorem 1.2 implies 

Theorem 1.3. Consider the points E^^\ I e Z, defined by 
(1.30) (/.i(E(')) =7r/2 + 7rl, 

and situated in Jg. Let L be the set of integers I such that E^''^ exists. Then, for 
sufficiently small e, 

• the distances between these points satisfy the inequalities 

ciS < \E^^^ - E^^-^^ \ < C2e, 

where Ci and C2 are two positive constant independent of e. 

• there exists a collection of intervals such that Ii C J 5 and such that 
the spectrum of equation (1-24) in J 5 is contained in Ui^lIi. 

Moreover for I G L, we have 

• the interval Ii is in an o{s)-vicinity of E^''\ 

• the length of Ii is equals to {t (E^^^) + ti {E^''^) (1+0(1)) . 

This theorem is similar to a result on the spectrum of Harper equation first 
proved in [14]. We prove it by using the idea on the relation between the Bloch 
solutions of (1.1) and (1.9) discussed above. 

1.4.2. Consider the monodromy matrix described in Theorem 1.2 for E being in 
one of the subintervals Ii described in Theorem 1.3. Clearly, the leading terms of 
its asymptotics have the form 



(1.31) M 
where 




±1 

Fi — 2X1 cos z 



Xi = t, (EW)A(EW), Fi = cos ME) /t{E^'^). 
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If we forget about the error terms and let x = I'Kipje + C, then equation (1.9) 
becomes equivalent to the Harper equation 

(1.32) / (x + /O + / (:r - h) - 2Xi cosxf (x) = Fi{E) f (x). 

This leads to a list of conjectures on the spectral properties of family of equations 
(1.1). Before discussing these conjectures, we have to make a remark about the 
individual equations of this family. 

1.4.3. For a general metrically transitive family of equations, the spectral prop- 
erties have to be formulated in probabilistic terms (see, for example, [20]). In 
particular, certain basic spectral results appear to be valid with probability one 
(i.e. for almost all the individual equations). Here, it is the parameter (/) that 
plays the role of the "random variable". Trying to study the spectral properties 
of one of the equations (1.1), i.e. the equation corresponding to a given value of 
e M, one has to restrict the solutions of (1.9) on the lattice 

Xn = nh — (/), n ^Ia. 

The restrictions of the solutions 

Xn = x{nh-(j)), neZ, 

satisfy the finite difference equation 

(1.33) Xn+i = M {nh - (f)) Xn, neZ. 

So, trying to analyze typical properties of individual equations of the family (1.1), 
one has to study solutions of equation (1.33). 

1.4.4. Consider the finite difference equation 

(1.34) tl^±l±^ + Xcos{nh + ip)^ljn = E4;^, 4' ^ h- 

It is the well known and intensively studied Almost Mathieu equation (for an 
almost complete list of references see [16], [15]). 

Clearly, (1.34) coincides with the restriction of (1.32) to a lattice. It is believed 
(and, actually, proved in many situations) that 

• if A > 1 and h is Diophantine, then the spectrum of (1.34) is pure point for 
almost all G M; 

• if < A < 1 and h ^ Q, then the spectrum is absolutely continuous for 
almost all G M. 

• if h ^ Q and A > 0, then the spectrum is a Cantor set. 

This immediately leads to three conjectures on the spectral properties of (1.1). 
However, in [13], it is proved that, for A = 1 and a full measure set of h, the 
spectrum of (1.34) is singularly continuous. But since the spectrum of (1.34) 
crucially depends on the value of A, the analysis of the spectrum of (1.1) on 
the intervals where A^ = ti{Ei)/t {Ei) is exponentially close to 1, or where the 
difference 1 — ti{E)/t{E) even changes its sign, has to be very delicate and, so, 
the results known for the Almost Mathieu equation lead here not to a conjecture, 
but to an open question. 
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Here, we collect some well known facts (see, for example, [17], [18]) on one- 
dimensional periodic Schrodinger operators (1-14). We assume that V is suffi- 



2. Periodic Schrodinger operators 

Here, we collect some well known facts (see, for example 
dimensional periodic Schrodinger operators (1-14). We assu 
ciently regular real valued function, 

V (x + l) = V (x), xeR. 
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The results described below are valid iiVE Lf^^{W). 

2.1. Bloch solutions. Let be a solution of (1.14) satisfying the relation 

(2.1) ^{x + l) = Xi){x), Va;eM, 

with A independent of x. Such solution is called a Bloch solution, and the number 

A is called the Floquet multiplier. 

Consider the spectrum of equation (1.14) in L^(M). It consists of intervals [Ei, E2], 
[E^, E4] . . . [E2n+i, £'211+2], . . . , of the real axis such that 

(2.2) El < E2 < Es < E4 . . . E2n < E2n+1 < E2n+2 < ■ ■ ■ 

(2.3) En -> +00, n -> +00. 

The intervals [i?2n-i5 -E'2n] {n > 1) are called the spectral zones, and the intervals 
{E2n, E2n+i), n > 1, are called the spectral gaps. Usually, the number of the gaps 
is infinite. 

Consider two copies of the complex plane E & C cut along the spectral zones. 
Paste them together to get a Riemann surface with square root branching points. 
We call his Riemann surface T. Denote one of the copies of this surface by r+ and 
the other one by r_. r+ (resp. r_) will be refered to as the upper (resp. lower) 
sheet. 

One can construct a Bloch solution 'tp {x, E) of the equation (1.14) meromorphic 
on this Riemann surface. The poles of this solution are located in the spectral 
gaps. More precisely, each spectral gap contains precisely one of these poles. All 
the poles are simple. They are located either on r_|_ or on r_. The position of a 
pole is independent of x. The Bloch solution ip [x, E) has the asymptotics 

(2.4) ^{x,E)=e'^''{l + o{l)), E^oo, 

(2.5) Eer, 0<arg£<27r. 

We define the branch of the square root so that, for E < 0, i\rE > if E' e r_|_. 

Denote by -0-1- the restrictions of the above defined Bloch solution on F-t. Outside 
the ends of the spectral zones, these functions are linearly independent solutions 
of (1.14). Along the cuts, we have 

(2.6) ^>bdx, E-iO)= ij^{x, E + iO), E21+1 <E< E21+2, I e N. 

2.2. Bloch quasi-momentum. 

2.2.1. Consider the Bloch solution tp{x, E) introduced in the previous subsec- 
tion. The corresponding Floquet multiplier A (E) is analytic on F. Represent it 
in the form 

(2.7) X{E) =exp{ik{E)). 

The function k (E) is called Bloch quasi-momentum. It has the same branching 
points as ip {x, E), but the corresponding Riemann surface is more complicated. 
To describe the main properties of k, consider the complex plane cut along the 
real line from Ei to +00. Denote the cut plane by Q. One can fix there a single 
valued branch of the quasi-momentum by the condition 

(2.8) iko{E) < 0, E < Ei. 
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The image of Co by ko is located in the upper half of the complex plane, 
(2.9) lmko{E)>0, EeCo. 

In Fig. 1, we drew several curves in Co and their images under the transformation 
E ^ kQ{E). 
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Figure 1. 

The quasi-momentum ko{E) is real along the spectral zones, and, along the spec- 
tral gaps, its real part is constant; in particular, we have 

(2.10) ko{Ei) = 0, ko{E2i±iO) = ko{E2i+i±iO) = ±7rl, I e N. 

Let Em be one of the branching points of k. Then both of the functions 

(2.11) f^{E) = {ko{E ± lO) - k^{E^ ± zO)) I^E-Era. E G M, 
can be analytically continued in a small vicinity of the branching point E^. 
Finally note that 

(2.12) A;o(£^) = v^ + 0(l/v^), |^| ^ 00, 

where E' G F and < avgE < 27r. This representation is uniform in argE' if 
< S < arg E' < 27r — 5 for any fixed 5. 

2.2.2. The values of the quasi-momentum ko on the different sides of the cut 
El < E < +00 are related to each other by the formula 



(2.13) 



ko{E + iO) = -ko{E - iO), Ei<E. 



Consider the spectral gap E21 < E < E21+1, I G N. Let Q be the complex plane 
cut from —00 to E21 and from E21+1 to +00. Denote by ki the branch of the 
quasimomentum analytic on Q and coinciding with ko for ImE > 0. One has 



(2.14) 



ki{E + iO) + ki{E- iO) = 2ttI, ii E < E21 ov E > E21+1. 



This relation follows from (2.10) and from the analyticity of the functions given 
by (2.11). 
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2.3. Periodic components of the Bloch solution. Fix I e NU {0}. Consider 
the solutions il)± introduced in subsection 3.1. They can be represented in the 
form 

(2.15) ^P^{x,E) = e'^'^l(^)''p^{x,E), EeQ, 
where p^{x, E) are functions periodic in x, 

(2.16) pf{x^l,E)=pf{x,E), VxeM. 

2.4. Useful formulae. We shall use the following two integral relations. 
Lemma 2.1. Let E ^ Ei, / e N, and let x e M. Then 

(2.17) / 'il)+{x,E)^_{x,E)dx = -%'k!{E)w{E), 

Jx-l 

where w = w {'ip+, ip-) and 

2 r ^{x,E)i;4x,E)dx = 

= -{ik"{E) + {k\E)f ) w (E) - 2ik\E) w (x, E), ^_(x, E) J . 

Here w (/ {x), g{x)) = f'{x) g {x) - g'{x) f (x). 

Proof. Though these two formulae are known, we describe their simple proofs 
here. Begin with the formula (2.17). Derivating (1.14) in E, we get: 

Multiplying this equality by ■0_ (x, E), integrating from x — 1 to x, then integrating 
twice by parts, we obtain 

(2.20) -w (^^{x',E),ij_{x',E)^ = ^Mx',E)ij4x',E)dx'. 

Now, by means of the formulae 

V^±(x-l,E) = e^^^ (^) Mx,E), 

we can express all the functions in the left hand side of (2.20) in terms of functions 
of the argument x' = x. This leads to 



w 



^(x',E),^_(x',E) 



= ik'{E)w{E), 

x'—x—l 



which implies (2.17). The formula (2.18) can be obtained similarly. But now, 
instead of (2.19), we start with the relation 



x.E) + V (x) -^{x,E) = E—^ (x, ^) + 2 — f (x, E) . 



9x2 5^2 V ' / ■ ^ ' dE^ ^ ' ' 9^2 V ' / ■ QE 

We omit the elementary details. □ 
In section 4, we shall also need the 
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Corollary 2.1. Under the assumptions of Lemma 2.1, 
(2.21) 

-w{^+{x,E), ^{x,E)^+J^ ^^+{x',E)^_{x',E)x'dx' = 

= -w'iE) + w (E) (-1 A In k'iE) + k'iE) - , 

where 

uj+iE) = -J^ ^ix,E)p_ix,E)dx Q p+ix,E)p_ix,E)dx^ . 

This corollary easily follows from formulae (2.17) and (2.18). Note that, 
though the Bloch quasi-momentum k is defined up to an integral multiple of 27r 
(and, therefore, the periodic components p± of 0± are defined up to some factor 
g=F27nix-j^ all the objects in (2.21) are uniquely defined. 



3. Main theorem of the complex WKB method 

In this section, we consider the family of equations (1.19), define the corre- 
sponding canonical domains and, then, prove Theorem 1.1. 

3.1. Canonical domains. The canonical domain is the main geometric notion 
of the complex WKB method. Its origin is tightly connected with the construction 
of the analytic solutions having standard asymptotic behavior in the parameter 
^ e C. 

3.1.1. The canonical domains can be described in terms of the complex momen- 
tum k((^). Remind that it is defined by formula (1.15) in terms of the Bloch 

quasi-momentum k iE) of the operator Hq = — — + V (u) with the same V 

du^ 

as in (1.19). The complex momentum k is a multi-valued analytic function. Its 
branching points are related to the branching points of the quasi-momentum by 
the relations 

(3.1) Ei = E -W i^p), leN, 

where Ei are the ends of the spectral zones of the operator Hq. The conformal 
properties of k and the structure of its Riemann surface depend on the analytic 
properties of the function W and on the value of the spectral parameter E. For 
W ((/?) = cos the analytic properties of the complex momentum and of its 
Riemann surface are described in subsection 4.1. 



3.1.2. In what follows, Ci and C2 ai'e two points in C such that 

Im < Im (^2 ■ 

These points can either be situated in the complex plane or be the points ±?cx). 
We shall denote by 7 a smooth curve going from Qi to C2- When we will integrate 
along such a curve, the curve will always be oriented from Qi to C2- We now state 
the definitions. 
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3.1.3. We call a smooth curve 7 C C vertical'ii it intersects the lines Im 2; = Const 
at non-zero angles 6. 

If there is a positive number such that, at any point of 7, the intersection angle 
6 satisfies the inequality 

(3.2) 5i<e <t:-5i 
then, we call 7 strictly vertical. 

3.1.4. Let 7 be a vertical curve. On 7, fix a continuous branch of the momentum 
K. We call 7 canonical if, along 7, Im J*^ k d(p is monotonically increasing, and 
Im /'^(k — tt) dip is monotonically decreasing. 

Assume that 7 is strictly vertical. If there exist two positive numbers 62 and S3 
such that 

(3.3) Im / Kd(fi > 62110. (ip' — (fi) 

J (p 

(3.4) Im / {k — t:) dip < —S^lm. {ip' — (f), 

J if 

then, we call 7 strictly canonical. 

3.1.5. Let K a simply connected domain containing no branching points of the 
complex momentum k. On K, fix a continuous branch k. We call K a canonical 
domain if it is the union of canonical curves connecting some two points Ci and 
C,2 fixed as above. To be precise, we have to exclude these two points from K. 

If there are three positive numbers di . 82 and ^3 such that if is a union of canonical 
curves satisfying (3.2) - (3.4), then, we call K strictly canonical. 

3.1.6. Assume that K C C is a canonical domain. Denote by dK its boundary. 
Fix a positive number 5^. We call the domain 

(3.5) C = {zeK\dist{z,dK)>5^} 
an admissible sub-domain of K. 

Note that the branching points of the complex momentum are outside of C, at a 
distance greater than 5^. 

Remark 3.1. We will describe the other classical geometric objects and construc- 
tions usually related to the complex WKB method (for example, Stokes lines) only 
in the case when M (93) = cos((^). These constructions will not be needed neither 
to formulate nor to prove our main theorem. However they are essential to un- 
derstand the geometry of canonical domains for a given equation. 

3.2. Canonical Bloch solutions. To describe the asymptotic formulae of of the 
complex WKB method for equation (1.1), one needs specially normalized Bloch 
solutions of the equation (1.21). 
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3.2.1. Assume that D is a simply connected domain in the plane cp G C contain- 
ing no branching points of k. 

Consider (1.21) as an equation with spectral parameter £ and consider the asso- 
ciated Riemann surface T. Let -0 {u, E) be the Bloch solution of equation (1.21) 
meromorphic on F. 

The mapping ip ^ E — W {(p) maps D onto a domain D C C Clearly. D is simply 
connected and does not contain any branching point of iJj {u, £). So, since F is a 
two sheeted Riemann surface, we can also consider two different copies Z)+ and 
l)_ of the domain D. 

Denote the restrictions of i> {u, £) on D± by ip±. These functions are two Bloch 
solutions of (1.21). Fix on a continuous branch of the Bloch quasi-momentum 
k{£). Clearly, 

^±(m, £) = e±^'=(^)«p±(M, 8), 

where p± are 1-periodic in li, 

p±{u + 1, £:) = p±(u, £). 

So, ±.k{£) are the quasi-momenta of '4>±{u, £). Note also, that p±{u, 8) are 
meromorphic m 8 E D. 

3.2.2. Let 

(3.6) -Ji£i^iitei|;^, s,o. 

Jo p+[x, 8)p-{x, 8)dx 

In view of (2.17), g± are meromorphic in D. 
Define 

(3.7) /±(«,£) = y^e/'^±(^)^^^° (^,£), 8^D. 
We call the functions 

(3.8) ■^±{u,v) = U{u,8), 8 = E-W{p>), 

canonical Bloch solutions associated to the domain D. Note that they are inde- 
pendent of the choice of the branch k. One has 

Lemma 3.1. The Bloch solutions ^'±(m, (p) are analytic in (p ^ D. 

Proof. It suffices to prove that the functions f±{u,8) are analytic in D. Let us 
prove it for /+. The case of /_ can be dealt with in the same way. In what follows, 
we denote by ~ the equalities up to some terms analytic in 8 & D. 

Obviously, the formulae from Lemma 2.1 remain valid for f±{u,8) when 8 & D. 
They imply that 

(3.9) S^i^)~-\i^'^m-^^J^^^ [?^±(u,S),Mu,S)) 



We have omitted the term | k'{8) since it is analytic outside the branching points. 
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Consider the second term in the (3.9). Remember that the poles of 'il)±{u,£) are 
independent of u. Let ^ -D be a pole of tp±{u, £). In a neighborhood of £q, we 
have 

where v± is the multiplicity of £o as a pole of ip±{u, S), and ip±{u, £) is a solution 
of (1.21) analytic in a neighborhood of £q. 

tp±{u,£) is a Bloch solution with the same Bloch quasi-momentum ±k {£) as 
ip±{u,£). Since k {£) G ttZ only at the branching points, the solutions 'tp±{u,£) 
are linearly independent. These arguments imply that, in a vicinity of £o: 



V 9£ 

and, thus, (3.9) takes the form 



u=l ^ - ^0 



2d£ ' ' £-£o 
This implies the analyticity of /-t in a vicinity of any point £q & D. □ 
In the sequel, we shall use the notation 

(3.10) a;±(^) = -W'{ip) ■ g±{E - W (^)). 

3.3. Formulation of the main theorem of the complex WKB method. 

We prove 

Theorem 3.1. Fix a positive U. Let K be a bounded strictly canonical domain 
for (1.19). 

For sufficiently small positive e, there exists a consistent basis f±{u,(p) for 
the family of equations (1.19) with the following properties: 

• for any fi/xed u G M, the functions f±{u, (p) are analytic in (p ^ C; 

• the functions f±{u,(p) have the asymptotic representations 

(3.11) f^{u,^) = e^s Io^d'P{-qf^{u,ip) + o{l)), e^O, 

-U <u<U, G C, 

where C is a fixed admissible sub-domain for K and ^± are the canonical 
Bloch solutions associated to the domain K. 

• the error estimates in (3.11) are uniform in (p and u. 

Moreover the asymptotic representations (3.11) are once differentiable in u. 

The rest of this section is devoted to the proof of this theorem. For the sake of 
brevity, we do not check explicitly the uniformity of errors nor the differentiability 
of the asymptotics: this does not require any additional estimates or discussions, 
and the reader will easily fill this gap. 



3.4. Sketch of the proof. Before beginning the actual proof of Theorem 3.1, let 
us outline the main ideas of this proof. 
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3.4.1. Let us assume that ip belongs to a bounded domain D C C and that u is 
in a fixed interval [—U, U] C M. Clearly, 



W {eu + if ) = W {ip) + 0{£), £^0. 

So, we may "replace" (1.19) by equation (1.21) with £ given by (1.22). Assuming 
that the domain D does not contain any branching points of the complex momen- 
tum K {(p) = k{E — W (<^)), we introduce two linearly independent Bloch solutions 
tp^{u,(p) of (1.21)-(1.22). Assuming that these two solutions have no singulari- 
ties in D, we prove that the initial equation (1.19) has two linearly independent 
solutions -ipl. close to -0° when e — >• 0. 

3.4.2. Though the constructed solutions tp]_{u,(p) are linearly independent and 
analytic in ip E D, they need not be consistent. If there exists any solution / {u, Lp) 
of (1.19) satisfying the consistency condition (1.20), it certainly can be represented 
in the form 

f{u, p) = A{ip) ijliu, p) + B{ip) ij^_{u, ip) 

where the coefficients A and B are independent of u. The consistency condition 
is equivalent to a difference equation of the form 

The coefficient of the matrix T can be calculated in terms of the Wronskians of 
the solutions tpl_{u — l,p + e) and tp}i-{u, p). 

3.4.3. Then, we compute the asymptotics of the matrix T to find out that 

■ ^iK(v)+0(s) 0{£) \ 

Using a transformation of the form 

A, (^) = e^i:,-i^)'i^+oie) ^^^^^ ^ j.^«(^)rf^+o(.) ^^^^^ 

we come to the equation 

(3.13) 

where the norm of the matrix Ti is small. 

3.4.4. The operator in the left-hand side of (3.13) can be easily inverted (see, for 
example [7]); so that, finally, we come to the equation 

where K is a singular integral operator acting in a space of Holder functions 
defined on a vertical curve j G D. If the curve 7 is strictly canonical then the 
norm of the operator is small. So, we can construct a solution of (3.14). This 
solution can be analytically continued in the domain D, and, as a result, one gets 
an analytic solution satisfying the condition (1.20). 
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3.5. Solutions of equation (1.19) close to Bloch solutions of equation 

(1.21). Let D be a bounded domain in the complex plane (p not containing any 
branching points of the complex momentum. Define, as in section 3.2.1, the 
functions ip±(u,8), and let 

(3.15) il^l{x,ip) = il^^{x,E-W{ip)), ip&D. 

In the sequel, we shall assume that D does not contain any poles of '0±(^) '^Y- 

The functions i>%{x, ip) are two linearly independent Bloch solutions of the equa- 
tion (1.21) (with S given by (1-22)) with the Bloch quasi-momenta ±k {(p) related 
to the quasi-momenta of ip±{u, £) by the formula 

k{^) = k{E -W {^)). 

We denote the Wronskian of the solutions "0° (x, (p) by wq[p)) i.e. 

= w {il)l{.,'p),il)%,ip)). 

Clearly, wq{ip) is analytic and has no zeros in D. 

Fix [/ > 0. Let us construct two solutions of equation (1.19) close to -0° for 
£ ^ uniformly m p E D and u G [— [/, U]. 

3.5.1. Fix a positive Ui so that < C/i < C/. Define an integral operator acting 
on C ([—[/, [/]) by the formula 

(3.16) 

U fU 





Rf{u,(p) = i^+{u,p>) / i)_{x,ip) f {x)dx -i)_{u,p>) / i)^{x,(p) f {x)dx. 

J-Ui JUi 

This operator coincides with a Green's function of equation (1.21) - (1.22). Indeed, 
letting 

Ho = -^ + V{u) + W{^)-E, 

we get 

HoRf = -wo f. 

Clearly, under our assumptions on V and W and the domain D, the operator R 
is bounded, 

(3.17) sup ||i?|| < C. 

3.5.2. For tp]., take the ansatz 

(3.18) V4 = ^± + £i?^± 

Substituting this ansatz into (1.19) and assuming that tp]. are solutions of this 
equation, we get 

(3.19) -A-i;i = -wog± + eA-Rg±, 



^From the very beginning, we could have considered the canonical Bloch solutions associated 
to the domain D (since they are analytic in D). However, we did not want to introduce these 
objects as an "ansatz" but show how they arise naturally in the constructions of the complex 
WKB method. This wih neither comphcate nor enlarge our analysis. 
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where 

(3.20) A.A(«.ri = II(f^i±^l^iIM. 

£ 

3.5.3. Since the Wronskian wq {(p) has no zeros in D, then 

, A , 

— < C, -U <u<U, (pe D, 

Wo 

where C is clearly independent of £, u and (p. Therefore, for sufficiently small £, 
equation (3.19) has a unique solution in C{[—U,U]); this solution satisfies 

(3.21) sup||^±||oo < C. 
The formula (3.19) implies that 

(3.22) g± = — il)l + 0{£), <peD, -U < u < U 

Wq 

in the sense of the supremum norm. 

3.5.4. By formula (3.18), we reconstruct the function in terms of ^'-t. This 
leads to the following results: 

1. for any fixed (p G D, tpl. satisfies (1.19) in C{[—U, U]); 

2. for any fixed u & [—U,U], is analytic in G 

3. if —U < u < U and ip & D, then 

(3.23) i;i{u,^)=i;l{u,^) + 0{e), 
the error estimate being uniform in (p and u. 

3.5.5. We shall need the asymptotics for the Wronskian wi (ip) of the solutions 
ip\{x, (p) and ip^{x, (p). Therefore we prove the 

Lemma 3.2. For £ — > 0, 

-u 



(3.24) w,{<p) = Wo{ip)+£W'{<p) I ' xPl{x,<p)xP^{x,<p)xdx + 0{£'') 

J-Ui 

uniformly in ip ^ D. 
Proof. By (3.18), 

= ^o + e {w{il^l,Rg.) + w{Rg+,il^'l))+0{e'). 
Now, using the definition (3.16) of i?, we obtain 

w {^iPj^, Rg ) + w [Rgj^, ^^_) = —Wq I ip+9-dx + Wo / ip^g^dx. 

Jui J-Ui 

Substituting the expressions given by (3.22) for g^. in this formula, we get 

w{i)l, Rg^) +w{Rg+, = / A i)t dx + O (s) . 

J-Ui 

This formula and formula (3.25) imply the desired result as 

A {ex, (f) = W'{(p)x + 0{£). 
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□ 

Note that, since wq has no zeros in D, Lemma 3.2 implies in particular that, for 
sufficiently small e, the newly constructed solutions are linearly independent for 
any ip ^ D. 

3.6. Consistency condition. The functions ip\.{u,ip) satisfy (1.19), but they 
need not form a consistent basis of solutions for the family of equations (1.19). 
However, in terms of these two solutions, we will construct two consistent ones. 

Let A[(p) and B[ip) be analytic functions of ^p. Define 

(3.26) 'il){u,^) = A{^)il)\{u,^)^B{^)il)\{u,^), <peD, 

Clearly, Vj(u,cp) verifies (1.19). This solution satisfies the consistency condition 
(1.20) if and only if, for Lp, cp + e G D, u eM., A and B satisfy 

(3.27) A{(p + £)i;l{u-l, ip + s) + B{ip + s)i;l{u-l, (p + e) 

= A (^) i^liu, ip) + B (ip) i;^_{u, ip). 

This relation implies that 

(3.28) (^i:::i)=rM(^l:;). 

where T is a 2 x 2-matrix with coefficients 
(3.29) 

_ w{xPl{.,^),iji{. -1,^ + e)) _ wixpl{.,ip),iji{. -l,<p + 8)) 

-ill — 7 ; ^ ) -'12 — 7 — ] — 7 , 

wi((p + e) wi(ip + ej 

_ w{i;l{. -l,^ + e),^l{.,ip)) _ w{i;l{. + ^i(.,y)) 

-'21 — 7 ; 7 , -1-22 — 7 ] 7 ■ 

wi((p + e) wi[(p + ej 

Now, our aim is to construct two analytic solutions of the difference equation 
(3.28) to get a consistent basis for the family of equations (1.19). Therefore, we 
shall first study the matrix T. 

3.7. Asymptotics for the coefficients of the matrix T. When constructing 
the solutions of the difference equation (3.28), we shall need to know the diagonal 
elements of the matrix T up to an error of size ^(e^), and the anti-diagonal 
elements up to an error of size O (e). 

One has 

Proposition 3.1. For sufficiently small e, 

Tn = exp (^iK {ip + e/2) + ^ ^ + ^^li'P)^ + O (e^), 

T22 = exp (^-iK {<p + e/2) + | ^ + ^^-i^P)^ + O (e^), 



ri2 = C»(e), T2i = 0{s), 
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where the error estimates are uniform in ip E D. Here, the functions uj^ are 
defined by formulae (3.10), and = k'{E - W {(pj) = (^J^k^ - W{ipj). 

Proof. The estimates of T12 and T21 are elementary. For example, as tp]. are 
analytic in e D, 

(3.30) 

T12 = 7 ^ = + O {£), 

and since 

the leading term in (3.30) is zero. 

To derive the asymptotics of Tn and T22, we shall use Lemma 2.1. We shall check 
only the formula for Tu. The second one is obtained in a similar way. 

In view of (3.18), 

^p\{u — l,(p + e) = ^p^{u — 1, (p + e) + e (Rg-) {u — l,(p + e) 

= e^^ ('^ + ^)^° («, ^ + e)+e{R g.) {u-l,ip + e) 
= e"^ («, ^) + s («' (cp) e"^ («, <p) + 

(^)^ ^) + ^Rg_) (u- 1,^)^+0 is'). 

Now, again using (3.18) for ijj^ and the result of the above calculation, we get 

w{-tpl{. ,(p),^l{. -l,(p + £)) = e'^'^Wo + £ (^iK' e'^'^ Wo+ 

(3-31) +e''^w{i;l{u), ^ (u)) + w {Rg+{u), e''^i;''_{u)) + 

o(p 



-w 



{^l{u),Rg.{u-l)))-rO{e'). 



Here, for the sake of brevity, we do not write the argument Lp anymore. Elementary 
calculations show that 



/u 
il)^_{x) g+{x) dx, 
■Ui 

w{iPl{u),Rg_{u-l)) = -Woe"^ / i^lix) g_{x) dx, 



'Ui 

and, using formula (3.22) as when proving Lemma 3.2, and the analyticity of W, 
we get 

w{Rg+{u),e'^r-{y'))+^{i^W, R9-{u-l)) = 

= e^>^W' ( f \l{x) ij'^_{x) X dx + f ijl{x)ij'^_{x)xdx] +0{e). 

\J-Ui Ju-l ) 
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Substituting this expression into (3.31) and taking into account Lemma 3.2, we 
obtain 

(3.32) 

+w{ijl{u),^{u)) + W' f ijl{x)ij'^_{x)xdx) +0{e^). 

(y'P Ju-l J 

But, in view of Corollary 2.1, we have 

wiil^liu), ^{u)) + W' tjjl{x)'4^°_{x)xdx = 

= u;o + ■iwo Q ^ In /cb - ^ k'(^) + u)%^ . 
Substituting this result in (3.32), one easily gets 

= e^'^ (^Wi+£ + w'q + Wq Q ^ In + t<j° j ) ) + ^ ^^^^ 

= e^'^ wi{ip + e) (^1 + £ (^k'I2 + 1 A in A;^ + ^ + O (e^) 

= wi{ip + e) exp (^K [if + e/2) + | ^ + ^^^4-) + ^ (^^)- 

We have indicated explicitly the cases where the argument of the functions differs 
from ip. This, together with the definition for Tn given by (3.29), implies the 
desired result. □ 



3.8. Transformation of the difference equation into an integral one. To 

construct analytic solutions of equation (3.28), we transform it into an integral 
equation. This transformation consists in three steps. 



3.8.1. Note that, in view of Proposition 3.1, the matrix T is close to a diagonal 
one. First, we transform the functions A and B so that the diagonal elements of 
the new matrix be close to 1. Therefore, we define 

(3.33) A((^) = e^^^(^)Ai(^), B{^) = e^^B{v) B,{^). 
So that, equation (3.28) takes the form 

(3.34) (^t::j)-n.)(^;j;:j), A 

with 




^-eA{ip + e) + 0a(v?) Tn e-^^i^ + ^) + t, 
^-eeiip + e) + eA{ip) T., e-^Bi<P + e) + 0s (v^) ^ 
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Fix a point (po in D and set 

(3.35) 9a{^) = - / Kdip + -\nkE+ co^dip 

(3.36) Ob{^) = — / Kdip + -\nkE+ uj^_d(p. 

^ J inn ^ J ipn 



In this case, in view of Proposition 3.1, 
(3.37) f = I + fi, Ti 



O (£2) o {e) 

2^0 O (e) O {s^) 



where 



(3.38) do = - I Kd(p, 

and the symbol O (e^) denotes functions satisfying the estimate 

\0{e%<Ce\ (peD, 

uniformly in if. 

3.8.2. A simple example of how to pass from a difference equation to an integral 
equation is given by the following simple 

Lemma 3.3. Let f (z) be a function of z e C analytic in the strip Fi < Im^; < Y2, 
Yi, Y2 G M. Fix Z\ and Z2 so that Imzi = Yi, lmz2 = Y2. Let 7^ be a smooth 
vertical curve going from z\ to Z2 via z. The function 

^/(^) = ^ / cot7r(C-^-0)/(C)c?C 
is a solution of the equation 

g{z + l)-g{z) = f{z) 
analytic in the strip Yi < Im z < Y2. 

The proof of this lemma immediately follows from the residue theorem. For anal- 
ogous statements see, for example, [8]. 

We will follow the idea of this lemma to transform equation (3.34) into an integral 
equation. 

3.8.3. We are ready to make the last two steps of the desired transformation. 
However, to construct the solutions /-t described in Theorem 3.1, one has to make 
these two steps in two different ways depending on whether one is dealing with /+ 
or /_. We shall concentrate on /+. The solution /_ can be constructed similarly 
and we shall discuss the details of this construction in the end of section 3. 

Fix two points Ci and C2 in D so that ImCi < ImC2 and so that they can be 
connected by a vertical curve 7 C D. For a function / {(p) sufficiently regular on 
7, define 

7r{(p' — (p — 0) 



(3.39) L+f{cp) = — / cot 



2i£ L 
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To construct a solution of (3.34) analytic in the domain 



(3.42) ' '"^1 ^ - ^ 1 




we consider the equation 

3.8.4. We need to make one more transformation. We fix a such that 

< a < 1, 

and let 

(3.41) 5l=£" 6^^0 52, 

so that equation (3.40) be transformed into 
where G 7 C D and 

(3.43) K^f (^) = {^) i^^29,{.) ^ 

In the next subsection, we shall show that, if 7 is a strictly canonical curve, then, 
after a proper choice of the functional space, the operator in the right hand side 
of (3.42) is small, and thus, that the equation (3.42) has a unique solution. Then 
we show that this solution can be analytically continued outside the contour 7. 
This will allow us to construct an analytic solution of the system (3.28). 

We note that an operator similar to was studied in [7]. It was the integral 
operator with the kernel 

acting in the space of functions analytic in a vicinity of a strictly vertical curve 
going from —zoo to +zoo. One can immediately see that, if, along this curve, 
Im P p dC, is strictly monotonically increasing and Im (p— tt) dC, is strictly mono- 
tonically decreasing, then \C (2, z')\ < Const ^-5\\uiz-\uiz'\ih ^ ^ > 0, and, thus, 
for small the corresponding operator is bounded. In our case, the estimates are 
more complicated since is a singular integral operator. 



3.9. Estimates of the norms of the integral operators. In this subsection 
we estimate the norms of L+ and as the operators acting in a proper functional 
space. We assume that the curve 7 defined above is sufficiently smooth. In the 
sequel, the letter C will denote a positive constants independent of e. 
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(3.45) ||/||a,6 = sup|/,(^)|+ sup 



3.9.1. Functional space. Let 7' be the curve 7 without the end points. Fix two 
positive numbers a and b so that 

(3.44) < 6 < a < 1. 

For a function / defined on 7', we let 

|/a(y)-/«(yOI 
I lib ' 

where 

/,(^) = iv.— cirk-C2r/(^). 

We say that / belongs to Ca,b{^]') if ||/||a,6 < co. One can easily see that 

1. The set Cab{l) equipped with the norm (3.45) is a Banach space. 

2. For any /eCi(7), 

(3.46) \\f\\a,b<C\\f\y 
and 

(3.47) l|/||a,.<C 11/11^ + C 11/11^-^ ll/'ll^. 

with some C > independent of /. 

We have chosen the space Ca,b{l) as, in particular, it is well known (see, for 
example, [19]) that the singular integral operator defined by the formula 

(3.48) Sf (^) = / -4^^^ V, ^ e 7', 
is a bounded operator in this space. 

3.9.2. Estimates of the norms . The operators L_|_ and are bounded operators 
acting in Ca,b{l)- The estimate of the norm of does not require any additional 
condition on 7. The estimate of the norm of is the principal point in the proof 
of Theorem 3.1. It is in this estimate where one has to use the definition of the 
canonical curves. First, we shall describe the estimates of the norms of these two 
operators, and then, we shall prove them. The estimates given below are rather 
rough, but we do not need any sharper ones. 

Begin with the operator L+. 

Lemma 3.4. Let 7 be sufficiently smooth strictly vertical curve. The norm of the 
operator in Ca,b{l) satisfies the estimate 

(3.49) ||L+|| < C£-^-". 

One can see that the norm of can be estimated from below be C / e. At the 
same time, if 7 is a strictly canonical curve, then the norm of satisfies a much 
better estimate. We recall that a and b are two arbitrary positive numbers fixed 
so that < 6 < a < 1. One has 

Proposition 3.2. Let 7 he suffi,ciently smooth strictly canonical curve. The norm 
of the operator in Ca^bil) satisfies the estimate 

(3.50) \\K+\\ KCe-"-''. 
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Proof of Lemma 3.4. For / e Ca,h(7), represent L+f in the form 

L+f = RLf + Sf, 
where S is the singular integral operator (3.48) and 



where 

1 / . . . . 1 



rL{<f) = - <f{ cot {tt(p) - i ) 

(p \ TT 

The norm of 5* is bounded by a constant independent of s. Let us estimate the 
norm of Rl. 

Note that, along any given strictly vertical curve passing via the point (po = 0, 

(3.51) \rLi<f)\<C, 
and 

f 1 \v\< 1, 

(3.52) K if)] < C i-27r \r]\ , ^ \n\ > 1 ^ = 1"^^- 

The constants here depend only on the parameter Si from the inequality (3.2) for 
a strictly vertical curve. Below, we shall denote Im ip by r]. 

First, we estimate \\Rl f\\c- Fix ip' G 7. Since the transformation cp does 
not change the angles on the complex plane C, the inequality (3.51) takes the 
form 



< C, (p, (p' e 7, 



e 

with C is independent of e, ip and ip'. This immediately implies that 
(3.53) ||i?,,/||c<C5-^ / ^ jr-J\f\\a,t = Ce-' 

Here 7^1,2 = ImCi,2- 
Similarly, using (3.51), we get 

\{RLfym<C8-' {I,{n) + h{v)+h{ri)) r; = Im^, 

where 

dr]' 



a,b- 



h 



T -^2 I H 

I2 — e 



|>e W - iv' - ViY iV2 - r]'Y 
m ^-2T,\r]-ri\lE 

\ri-n'\'>e iv' - viY im - v'Y ' 

Now, elementary estimates show that 

sup |/j|<Cei-", J = 1,2, 3. 

7)1<7)<7)2 
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Thus, 
(3.54) 



\\{RLfy\\c<Ce 



-l-al 



\a,b- 



In view of (3.46), (3.53) and (3.54) imply the desired result. □ 

Proof of Proposition 3.2. As in the proof of Lemma 3.4, we represent K+f in 
the form 



where 



K+f = RKf + Sf, 



rK {(p, (p') 



e {if, if') , 

£ TT 



, . A^'-^) .\ ^ r' K{s)ds 

e {(f, (f') = cot — - I \ ■'f ^ ^ 



We begin by estimating As before, let r] = lmf. Note that 



K (s) ds 



<C\(f-(f'\, (f,(f'e-f, 



Therefore, 

if-(f 

where the error estimate is uniform in e, f and f'. So, 



- (- + -0(^-^') K ^,^'^7, \r]-ri'\<e, 



(3.55) 



Vxiv, <C, \ri-'n'\<e, if.if' E'^f. 



Now, we assume that 7 is strictly canonical. Due to the definition of a strictly 
canonical curve, we show that 

(3.56) |e(^,^')l <Ce-f l^-^'l, - V| > e, y^, G 7, 

where 5 is a positive constant independent of e. 

First, we consider the case where r] — rj' > £. In this case, since 7 is strictly 
vertical, 



7r{f' - (f) . 
cot I 



< Ce~~ 



(r? - v') 



(3.57) 
and, therefore, 

I / /M -Im I"^,(k (s) — it) ds , 

In view of the definition of a strictly canonical curve, 

Im / (k (s) — tt) (is < —82 {rj — rj'), 
Jip' 

with some positive constant 82 (see (3.3)). This implies (3.56). 
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Now, we consider the case where r] — rj' < —e. In this case, instead of (3.57), we 
get 



cot I 



and, therefore, 

-- Im K is) ds 
But, since 7 is a strictly canonical curve, 



pip' 

Im / K (s) ds > S {rj' — rj) , 

J in 



with some positive (5, which again leads to (3.56). 

The estimate (3.56) and the definition of Tk imply the inequality 

5 u ^/i Ce 

|?7 — ' 

The inequalities (3.55), (3.58) and the definition oi Rk imply that 



(3.58) |rK(^, ^')l < + 



\r\-r\\>e, ^,^'€7. 



where 



\RK!{^)\<C{h{ri)-rh{ri)^h{r])) e 7, 

dri 



J m , I' 



1 



\<e iv' - Vi)" {V2 - n'Y 
drj' 



m, \v-v'\>s W - V\ iv' - ViY iV2 - ri'Y ' 



^-5\ri-r]'\/e 

£ Jm, \^-^'\>e - ViY iV2 - V'Y ' 
Elementary estimates show that 

sup |/J•|<Ce-^ j = 4,5,6. 

7)1<7)<7)2 

And, result, we get 

(3.59) \\RKf\\c<Ce-^\\f\U. 

Having estimated we estimate ||(i?ii-/)'||(:;. We first show that, for ^p, ip' e 

7, 

1 when \r] — rl\ < s, 



(3.60) 



iv, V') 



2^1^-^''/^ + ^, when|r?-V|>£, 



Note that 
(3.61) 



dip 
where 



i<^, ^') 



( ^ (^"1^ C (cot C-i)) ^ ~ 7 ^ v') 

_ 7r{p' - p) 



7r((^ — ip'Y ' 
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Now, first, we assume that I?? — r/| < e. We recall that 7 C -D, and, thus, the curve 
7 is situated at some distance from the branching points of the quasi-momentum 
K. Therefore, 

j K{s)ds={(p' -(p)K{(p) + 0{{<^-<^'f), 99, (^'€7, 
with the error estimate uniform in ip' and s. This allows us to get 

^ <-') = ^ + - ^ + ° (^) ) ■ 

(p,(p' e 7, |r/ - r/l < e. 

And, substituting this representation in (3.61), one easily comes to the first in- 
equality of (3.60). 

To treat the case where \ j] — rj'\ > e, one notes that, for It] — t]'\ > £ and (p,ip' e 7, 
one has 

I sin^C(cot C-i) 1"^ < C. 
This estimate and the estimate (3.56) lead to the second inequality from (3.60). 
Using the estimates (3.60), one obtains 

\{RKfn^)\<C6-' {h{v)+l2{r]) + h{v)) \\f\\a,b, ^€7, 
with the same /i, I2 and I3 as in the proof of Lemma 3.4. This gives the inequality 

(3.62) \\{RKfy\\c<Cs-'-\ 

similar to (3.54). Now, combining (3.59) and (3.62), by means of (3.47), we get 

\\RKf\\a,b <C{e-'^ + e-<'-'h-^'^'^^') \\f\\a,t < Ce-'^-''\\f\U,t. 
This implies the announced estimate of □ 

3.10. Analytic solution of equation (3.28). We are now ready to solve equa- 
tion (3.28). In this subsection, we prove 

Proposition 3.3. Let K be a strictly canonical domain related to some points 
Ci, C2 G C, ImCi < Im(^2- Let C he an admissible sub-domain of K. Assume that 
C does not contain any pole of ■0^ . For sufficiently small £ > 0, there exists a 
solution of (3.28) analytic in the strip 

(3.63) Im Ci < Im < Im C2. 
This solution has the asymptotics: 

Here, /5 is any fixed number, 

< /? < 1, 

and 9a is the function defined by (3.35) The error estimate is uniform in (p ^ C. 
The rest of this subsection will be devoted to the proof of Proposition 3.3. 
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3.10.1. Let the domain D in all the previous constructions coincide with K, and 
let 7 G K be one of the strictly canonical curves connecting and (2- First, 
construct a solution of (3.42) with the components Ai and B2 from the class 
Ca,b{j)- In view of Lemma 3.4 and Proposition 3.2, for sufficiently small e, the 
components of the matrix integral operator from the right hand side of equation 
(3.42) satisfy the estimates 

||L+ (O (£2) . )|| < Ce^-", ||L+ (O (£i+") . )|| < C^"-^ 

(3.65) 

\\K+ {O (e^-") . )|| < Cel-"-«-^ \\K+ {O (e^) . )|| < Ce^-a-fc^ 

where || . || denotes the operator norm induced by || . \\a,b- Fix the numbers a, a 
and b so that 

(3.66) 0<b<a<a<l, a + a + b<l. 

Clearly, in this case, for sufficiently small e, equation (3.42) has a unique solution 
in Ca^b- This solution satisfies 

yli = 1 + C»(e"-"), = C»(£i-"-"-''). 

Here the error terms are estimated in the || . ||a,6-norm. Letting 

(3.67) /5i = min{a — a, 1 — a — b} 
we get 

(3.68) Ai = l + 0{e^'), B2 = O {e^'-'') . 
Note that < /5i < 1. 

3.10.2. Now, let us check that Ai and B2 can be analytically continued in the 
strip defined by (3.63). Clearly, the poles in ip' of the kernel of the operator 
are situated at the points (p + + £/, / e Z. Therefore, the right hand side in 
(3.42) can be analytically continued in the strip 

S, = {ipeC: 3C e 7', < Re (y? - Re C < £} , 

where, as before, 7' is the curve 7 without the endpoints. Deformating the in- 
tegration contour in (3.39) inside S^, we see that, in fact, Ai and B2 can be 
analytically continued in the strip 52^. Repeating these arguments, we prove that 
these functions can be continued in the part of the strip defined by (3.63) situated 
at the right of 7. 

We check that Ai and B2 are analytic at the left of 7. Denote the second term in 

(3.42) by £ ^ ^ ^ [ip). We assume first that ip is located at the right of 7 and 

deform the integration contour so that it passes between the points ip and ip + h. 
Since the kernels of L+ and have simple poles at cp' = ip, instead of (3.42), 
one obtains 

(3.69) [/ + 0(e-)](^;(^;) = (J).£(^j)(.). 

Here, the second term in the brackets in the left hand side is a 2 x 2-matrix with 
analytic coefficients which can be estimated by Ce^~" for (p G K. This term 
is the residue that appeared in the course of the contour deformation. We have 
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continued Ai and B2 to the right of 7 by using (3.42). Applying almost the same 
arguments to (3.69), one continues these functions to the left of 7. 

3.10.3. Now, we discuss analytic properties of the functions Ai and B2 in some 
small vicinities of the end points of the curve 7. Clearly, for 99 G 7 being in a 
sufficiently small vicinity Vi of Ci, equation (3.42) can be rewritten in the form 

where S is the singular integral operator (3.48), and Ci is a matrix integral oper- 
ator with some kernel analytic in G Vi and in ip' G 7. Therefore, 

where / is analytic in Vi. We recall that Ai and B2 belong to Ca,b{l)- The analytic 
behavior of the singular integral in a vicinity of the ends of the integration contour 
is well known (see, e.g. [19]). Using classical results, we see that (i is a branching 
point of Ai and B2, and that these functions admit the representations 

(3.70) 

(^ - Ci)'^ A, (^) = ci + O (1^ - Cil'O, - CiT B2 iv) =C2 + 0{\^- 01"^, 

where ci^2 are some complex constants, and hi is a positive number. Similarly, in 
a sufficiently small vicinity V2 of C2, 

(3.71) 

(^ - C2)" Ai (^) = C3 + O (1^ - C2|'^), (^ - C2)" B2 (^) = C4 + O (1^ - (2!"^). 

We can assume that < 61 = 62 < ^- We set 

(3.72) /5 = min{a — a, 1 — a — bi}. 
Clearly, < /5 < 1. 

3.10.4. Consider the functions Ai and B2 along some other strictly canonical 
curve 7i G K. We show that, for sufficiently small e, they have the asymptotic 
representation 

(3.73) Ai = l + 0(e^), B2 = 0(e^-") in CaMili) 
similar to (3.68) 

Therefore, instead of (3.42), we consider the analogous equation along the curve 
7i. For sufficiently small e, this new equation has a unique solution in Ca,bi{"ii), 
and the new functions Ai and B2, e.g. the components of this solution, do have 
the desired asymptotics. But, the representations (3.70) - (3.71) for the "old" 
functions Ai and B2 imply that they belong to this functional space. And just 
deforming the curve 7 to 71, one can see these two "old" functions satisfy the 
"new" equation. This implies the desired result. 
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3.10.5. Consider an admissible domain C of K. Now, we can easily show that 
(3.74) Ai(^) = l + 0(e^), B2(¥^)=0(e^-"), ^ e C, 

with the uniform error estimates. 

Let (po ^ C. One can consider ipo as an internal point of a domain Kq <Z K 
with the boundary Fq consisting of two segments of two strictly canonical curves 
7i and 72 situated in K. We represent the functions Ai and B2 by the Cauchy 
integrals along Fq. Along Fq fl 71 (resp. Fq fl 72), the functions Ai and B2 admit 
the representation (3.74) in Ca,bi (71) (resp. Ca,bi (72))- This leads to the desired 
asymptotic formulae for any fixed (p E Kq. Moreover, the error estimates are 
uniform in 99 in a sufficiently small fixed vicinity of ipo. 

Thus, we see that any point (po of the compact C has an open vicinity Vq such that 
the desired asymptotic formulae are uniform m ip ^ Vq. This implies the needed 
result. 

3.10.6. To finish the proof, we reconstruct in terms of B2 the function Bi by 
formula (3.41). The vector ^ ^ {(p) satisfies equation (3.40). By means of the 
residue theorem, one checks that 

L^f, ( ^; ) (^+£) - L^f, ( ^1 ) (^) = ( ^; ) (^), 

Im <lunp < Im C2. 
This implies that this vector satisfies (3.34), and, by means of (3.33) we reconstruct 
a vector solution ( ^ J of the difference equation (3.28). Clearly, this solution is 



^ B ^ 

analytic in K; moreover the asymptotic formulae (3.74) and the relations (3.33), 
(3.41) imply (3.64). □ 

3.11. Consistent basis solutions. In this subsection, we complete the proof of 
Theorem 3.1. 

3.11.1. First, we construct the solution /+ described in Theorem 3.1. Denote 
the components of the vector solution of (3.28) described in Proposition 3.3 by 
.4+ and B^. Let tp\. be the solutions of (1.19) constructed in subsection 3.5. We 
consider the function 

(3.75) U{u, ^) = A+{^) i^liu, ^) + B+(^) i;l{u, ^), 

where ImCi < Im(^ < ImC2, m G M. This function is a solution of (1.19) analytic 
in (/? e C. By construction, it satisfies the condition (1.20) 
Finally, we note that the asymptotics (3.64) and (3.23) imply 

/+(«,^) = e^-(^) «(n,^) + 0(£^)), ^eC, -U<u<U. 

Remembering the definitions of 9a, "0° and and equation (3.35), we see that 
this formula coincides with (3.11). 
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3.11.2. Up to this moment, we were assuming that C does not contain any pole 
of -0° . Now, we get formula (3.11), and see that the leading term in this formula 
is the canonical Bloch solution. By Lemma 3.1, it has no poles. Having found 
the canonical Bloch solutions, we could have used them everywhere instead of 
This remark enables us to remove the assumption about the absence of the poles 
of ^pl in C. 

3.11.3. Let us discuss the construction of the second solution of (1.19) described 
in Theorem 3.1. It is similar to the one leading to /+. So, we shall only single out 
its new elements. 

Now, the central objects of the analysis are the operators 

L-f (^) = ^ / (cot ~^ ~ + / W) 



and 



acting on the "old" space Ca^bil)- As in subsection 3.9.2, one checks that, for 7 
being a strictly canonical curve. 

To reduce the analysis of (3.34) to the analysis of these two operators, we consider 
instead of (3.40) the equation 

^ \ I \ + L_T, ( ' ], ^ e 7 C K. 



Note that it differs from (3.40) not only by replacing L+ by L , but also by 



replacing the term ( n ) ( ? 1 ■ latter allows to obtain, by means of the 



transformation 



the integral equation 



for (/? G 7 C K. 

The analysis of this equation is similar to the analysis of (3.42), and, as a result, 
one constructs the solution /_ of (1.19) described in Theorem 3.1. Omitting the 
details, we note only that 

(3.76) /_(n, cp) = A-{cp) ijliu, ip) + B-{ip) ijl{u, cp), 
for Im (i <lunp < Im (2, u eM.. 

The functions A^ and B are analytic in the strip defined by (3.63) and admit 
the asymptotic representations 

(3.77) A-{^) = e^B iv) o (ef^), B-{^) = e^^ (v) (1 + (e'^)), 

for (p G C when e ^ 0, with the error estimate being uniform. Here, /5 is any 
fixed number, < /5 < 1, and Ob is the function defined by (3.36). □ 
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3.11.4. Finally, we check that the Wronskian of f±{u, ip) is independent of cp, 
and, thus, that they do form a consistent basis of solutions for the equation family 
(1.19). By (3.75) and (3.76), 



w 



{U,f.}={A+i^)B-{^)-A-{^)B+{^)) w,i^). 



Here wi{(p) is the Wronskian of the solutions '0^. Now, using the asymptotics of 
A'^ and described in Proposition 3.3, the asymptotics of A^ and B given by 

(3.77) , and the asymptotics of wi (see Lemma 3.2), we get 

(3.78) {/+,/_} = «Jo(^) exp (^\nk'E + J\io^ + iOo) difi^ (l + o(l))- 
By (3.6) 

^0 +^0 = - g^ In J -0° (x, if ) -0° (x, ip j) dx. 

Using Lemma 3.2, we transform this expression to the form 

d 

u)^ = -— In ( U'e WQ{ip) ) . 

This allows to get finally 

w {/+, /_} = U'e (1 + ^ (^)), g{^) = o (1). 

This formula shows only that the leading terms of the asymptotics of the Wron- 
skian is independent of (p. 

Let us discuss the term g{^p)) in the last formula. Since f± satisfy (1.20), it is 
e-periodic in ip. Moreover, its estimate is uniform in 99 e C. 

Let us change the notations, denoting by /+ the "old" function /+ divided by 
1 + g. The "new" functions /±, clearly, have all the properties announced in 
Theorem 3.1. This completes the proof of the theorem. □ 

4. Main constructions of the complex WKB method 
for equation (1.24) 

Here we consider the family of equations (1.19) with W {(p) = cos (p. 
4.1. Complex momentum. 

4.1.1. We begin by describing the Riemann surface F of the complex momentum. 
We recall that the conformal properties of k and the structure of its Riemann 
surface depend on the spectral parameter E and on the lengths of the spectral 

zones and the spectral gaps of the operator Ho = — — + V. For the sake of 

dx^ 

definiteness, we shall work under the condition (1.26) 

The branching points of the complex momentum k. are described in formula (3.1). 
As cos pi is periodic, each of the Ei gives rise to an infinite lattice of branching 
points of the complex momentum. There are two branching points corresponding 
to each of Ei and situated in the strip < Reip < 27r. Let us describe these 
branching points. 
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The branching points corresponding to Ei are real. They are symmetric with 
respect to the point tt. We denote them by (pi and 2% — (pi where < (pi < tt. 

The branching points corresponding to Ei, I > 1, are situated on the line Re(p = 
TT. They are also symmetric with respect to the point tt. The points with positive 
imaginary part will be denoted by (pi. Let 77^ = lm(pi > 0. 

The whole picture of the branching points is 27r-periodic. It is partly shown on 
Fig. 2. 

The conformal properties of the complex momentum near the branching points 
are determined by the behavior of the Bloch quasimomentum k near the points 
El (see subsection 2.1). In particular, the branching points ipi and 27r — (^i related 
to El are of square root type. 




4.1.2. One can fix a single valued branch of the complex momentum in a vicinity 
of (/? = by the conditions 



(4.1) 



ImKo(0)>0, ReKo(0)=0. 
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All the other branches defined there are then related to by the formulae: 
(4.2) i^mi^) = ±i^o{^) + 27rm, m e Z. 

The branch Kq can be analytically continued along any curve not passing through 
any branching point. In Fig. 3, we draw several curves on the complex plane 
(fi & C, and then, we show the images of these curves by Ko{(p). The results shown 
in these figures easily follow from the discussions of subsection 3.2. 

Note that on the shadowed domain in Fig. 3, we have 

Ko{(p) = ko{E - coscp), 

where ko is the branch of the Bloch quasi-momentum defined in section 2.2.1. 
To get the following four pictures, one uses known elementary properties of the 
conformal mappings. We shall use Fig. 4 for drawing the cuts on the complex 
plane needed to define the branch kq globally. 
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Cut the complex plane e C as in Fig. 2, e.g., first draw the horizontal cut 
between the points (pi and 2% — (pi, the vertical cuts between the pairs of the 
points ip2i and (p2i+i, I £ N, and then between all the pairs which can be obtained 
from the above ones by the reflection with respect to the real line and by the 
27r-translations. Denote the complex plane thus cut by C^q . 
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The branch can be continued up a single valued function defined on C^q . This 
functions has the following properties: 



(4.3) 
(4.4) 
(4.5) 



Ko{(p + 27r) = -Ko{(p) 



for e Qo . 

Relation (4.3) follows from the definition of Kq, formula (4.4) is a consequence of 
the observation that Ko{(p) G iR+ if < ip < (pi. The last relation comes from the 
analysis of Fig. 5. 
We finish this subsection with the 



Lemma 4.1. // 

(4.6) 

then, we have 
(4.7) 



Re ^ TT + 27rm, m e Z, 



Ko{ip) = e^^^/' + O (e-l^l/'), Im ^ ^ ±oo. 
v2 



This representation is uniform in Re (p outside any fixed 5-vicinity of the lines 
Re = TT + 27rm. 

Proof. The asymptotics (2.12) and the relation between ka^E) and Kq[lp) for 
being in the shadowed domain in Fig. 3 imply the asymptotics for IraLp ^ +oo 
in the case where < Ret/? < tt. Relations (4.3)-(4.5) show that this asymptotics 
remains valid for all Re<^ satisfying the condition (4.6). □ 
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4.2. The Stokes lines. The definition of the Stokes lines is fairly standard. The 
integral J Kd(p has the same branching points as the complex momentum. Let 
(po be one of them. Consider the curves beginning at (po and described by the 
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equation 

(4.8) Im r{K{0-K{^o))d^ = 0. 

These curves are the Stokes lines beginning at ipo. It follows from (4.2) that the 
Stokes line definition is independent of the choice of the branch of k in (4.8). 
Assume, as before, that E satisfies (1.26). In this case, at any branching point 
begin exactly three Stokes lines. The angle between any two of them at this 
point are equal to Some of the Stokes lines are finite: they connect pairs 
of finite branching points; some of the Stokes lines are infinite: they go from 
finite branching points to the infinity. The asymptotes of the infinite lines can be 
described by means of Lemma 4.1. 

The character of the Stokes lines in the domain < Re 2; < 27r is shown in Fig. 
6. 

Consider the branching point ipi. Two of the Stokes lines beginning at (pi go to 

±ioo and have the asymptote Reip = 0. The last Stokes line starting at (pi is 

finite; it goes along the real axis to 27r — (^1. In the same way, for each / > 1, two 

of the Stokes lines beginning at the point cpi are infinite. And the last one is finite 

and runs along the line Re<^ = tt from one branching point to another. 

The whole picture of the Stokes lines in the domain < Re 99 < 27r is symmetric 

with respect to the real line as well as with respect to the line Re 99 = tt. 

If is a branching point for k, then so is the point ip + 2ti. The Stokes lines 

passing through 99 + 27r are just the 27r-translates of the Stokes lines beginning ^p. 

4.3. Canonical domains. Here we describe some canonical domains in the case 
when W[lp) = cos{(p). 



4.3.1. In Fig. 7, we show two canonical domains. The first one (part A) corre- 
sponds to the case where (i = +ioo and (2 = —ioo- Part B corresponds to (i and 
(2 finite. 




consists of segments of Stokes lines. This domain exists if E' — 1 < Ei. 
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To check that Kq is canonical, one has to investigate the families of the lines 
Re Ko d(f = Const and Re /^(/«o — t^) dip = Const inside this domain. 

In Fig. 8, we show these families of lines. In the left hand side part of Fig. 8, we 
show some of the lines Re kq d^p = Const, and, in the right hand side, we show 
some of the lines Re /"^(^o — tt) dip = Const. The arrows indicate the directions 
of increase of the functions Im Kq dp and Im /'''(kq — vr) dp. 
In Fig. 9, we show one more canonical domain Ki. This domain exists if Ei < 
E + 1 < E2. It is canonical with respect to the branch fixed by the condition 
< Ki(7r) < TT. We note that Ki has the following symmetries: for p ^ Ki, 

(4.9) Ki{2tt-p>) = Ki{p), 

(4.10) Klip) = i^iip). 

The function ki can be continued on the complex plane cut along intervals [2k7r — 
Pi, 2k7r + pi] along the real line and as in Fig. 2 outside of the real line. In this 
case we have 

(4.11) Ki((^ + 27r) = -Ki(^). 




Relations (4.3) - (4.5) and (4.9) - (4.11) show that any domain obtained of Kq 
or Ki by means of the 27r-translations and the reflection with respect to the real 
line are canonical. 

4.4. Admissible sub-domains. We first note that the admissible sub-domains 
corresponding to Kq and Ki are bounded. Moreover we do not need the notion 
of strictly canonical domains anymore. Indeed, one can show that, if if is a 
canonical domain and A is one of its admissible sub-domains, then there exists B 
a strictly canonical domain such that A C B C K . Thus, the asymptotics given 
by Theorem 3.1 are valid on any fixed admissible sub-domain of a given canonical 
domain. 

41 




Figure 9. 




5. The computation of the monodromy matrix for equation (1.24) 

This section is devoted to the computation of the monodromy matrix asso- 
ciated to the consistent basis of solutions of (1.24) constructed in Theorem 3.1 for 
the canonical domain Kq defined in section 4. Before we begin, let us underline 
that we will commit a slight abuse of language and write the monodromy matrix 
M{(p) instead of M(</)) = M{(p/e). 

We now state the main theorem and then use the rest of the section to prove 
it. To describe the monodromy matrix, we define the following action integrals 

k \ , ti = exp ( 2- / (k — tt) 

(5.1) V - y V 



t = exp 2- 

£ .10 



K 



and 



2^ 



rV2 



Here the functions k, uj± denote the branches of these functions first defined in Kq 
and then continued along the path 7 (shown on Fig. 10) to the domain KqU Ki. 
We prove the 



Theorem 5.1. Define coefficients of the monodromy matrix M{(p) by 

Then, for any fixed S > 0, for e > small enough, uniformly for |Im(^| < 
Im (fis — S, we have 



M(^) 



t-e^<^i(l + o(l)). 



Jcj) 



ipJo Jo 



2z7r 



e'1'^{l + o{l))-ti-e e 



((y5 + TT + (i)2) 



(1 + «(!)) 



42 



'm2i{(p) = -ie 



/ + / 

JO ^0 



2m 



m22(^) = + 0(1)) + ^e^'^Hl + o{l)) 

Theorem 1.2 is an immediate corollary of Theorem 5.1 if one sets 

and if we prove that ^2 is real. This is done at the end of the proof of Theorem 5.1. 

For the canonical domain Kq (defined in the previous section), denote by (/^) the 
solutions of equation (1.24) constructed in Theorem 3.1. The monodromy matrix 
associated to this basis is then defined by 

F\u, ^ + 27r) = M(^)F°(m, ip) 



(5.2) 
where 



5.1. The basic ideas guiding the proof. In a subset of Kq where we know 
both the asymptotic behavior of F^{u,ip + 27r) and F^{u,ip), it will be easy to 
compute the asymptotics of M. Unfortunately a quick look at section 4.3 shows 
us that Kq n {Kq + 27r) = 0. Therefore we need to introduce an auxiliary canonical 
domain. 




5.1.1. Auxiliary canonical domains and transition matrices. The first auxiliary 
domain we need we have already defined: it is Ki (see Fig. 9 and subsection 4.3). 
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Let us denote by (/_{_) the consistent basis constructed by Theorem 3.1 in Ki. As 
(/£) and (/^) are both basis of solutions of (1.19), we define 

(5.3) F\u,if) = T^{ip)F\u,^). 

We will be able to compute the asymptotic behavior of Ti in Kq fl Ki . We will see 
that we will actually be able to go beyond this simple remark using analyticity, 
periodicity and a continuation lemma. 

The second auxiliary canonical domain needed will be denoted by K2. It is defined 
to be K2 = Kq + 27r. This domain is canonical for the branch of k obtained by 
continuing Kq along the curve 7 shown on Fig. 10. Moreover it is immediate that 
if we define f±{(p) = f±{^ — Stt), then for K2, this defines a consistent basis. As 
(/I) and (/^) are both basis of solutions of (1.19), we define 

(5.4) F\u,^) = T2{^)F\u,^). 

As above, we will be able to compute the asymptotic behavior of the matrix T2 
in KiH K2, and beyond that region. 

Putting (5.3) and (5.4) together, we get that 

F°(n, ^) = T,{^)T2{^)F\u, ^) = T,{^)T2{^)F\u, ^ - 27r). 
Hence the monodromy matrix reads 

(5.5) M(¥5) =ri((^ + 27r)r2((^ + 27r). 

5.2. The asymptotics of Ti and T2. These asymptotics are given by the 
Proposition 5.1. Let us define the coefficients of Ti by 



ci((^) fii(^) 



Then, for any fixed < F < Imips, for e > small enough, uniformly for 
\lm(p\ < Y, we have 

(5.6) ai(<^) = e-^So-°'i^+So<'i^{l + o{l)), 

(5.7) bi{^) = ei^o'''°'^^ -oil), 

e^(^) = -ie-Uo^od^ ^g2./;^.o<i^^j-.o,^_;.,^o,,^^ ^^^^^^ 

(5.8) ^ 

_e-27/r(«o-^)<iVgr2^°<iv-/r^V^e-2'^^(^-^)(l + 0(1))) , 

(5.9) di{^) = e-e^o"«od¥'+/o"-°^^(i + 0(1)). 

Here all the integrals are taken along contours living in Kq U Ki U K2. 
and by the 

Proposition 5.2. Let us define the coefficients ofT2 by 



C2i<f) d2i<f) 
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Then, for any < F < Ivmp'^, for e > small enough, uniformly for \lv[np\ < Y , 
we have 

f2T ^ J,„ I fSvi- ,.,1 . 



(5.11) 



(5.10) a2(^) = e^ ^'^""1^^+^'^" "+^^(1 + 0(1)), 

(5.12) C2(^) = e-^^^'""''"^'^^ -0(1), 

(5.13) (kiip) = e-eS'^''^'i^+S'^^-'i^{l + o{l)). 

Here all the integrals are taken along contours living in Kq{J KlV^ K2. 

The next subsections will be devoted to the proof of Propositions 5.1 and 5.2. 
The proofs being very similar we will not dwell on the one of Proposition 5.2, but 
just give the necessary modifications of the proof of Proposition 5.1 

5.3. The proof of Proposition 5.1. By definition, we have 

(5.14) a,{^) = ^w{flf'_), h,{^) = -^w{flJl), 

(5.15) ci(^) = -^(/°,/l), d,{^) = -—w{f'_Jl), 

(5.16) ..=.(flJi)=f.lfl-fl^Jl 

We recall here that (/±^) are the solution of (1.19) constructed in Theorem 3.1 for 
the canonical domains Ko,i. In admissible sub-domain of Ko,i: by definition, 
/_!_' have the following asymptotics when e ^ 

(5.17) /°'^(«,^) = e^^^'o.i'^M'i^(*°/(«,^) + o(l)) 
where ipo = 0, (pi = tt, 

and (resp. u]_) are the analytic continuations of -0° (resp. Wj.) along 7 from 
Ko to Ki (see Fig. 10 and 12). 

This implies that 

(5.18) w, = w{fl, fl) + 0(1) = k'{E + 1) ■ woi-ir) + o(l) 

where Wq is the analytic continuation of the Wronskian w{f^, /°) from Kq to Ki 
along 7 (see Fig. 10). 

5.3.1. Periodicity. One remarkable fact about {f±^) is that 

f'/{u,p + s) = f'/{u + l,^). 

Using the fact that the Wronskians taken in (5.14) are independent of the point 
u, we get that these Wronskians are e-periodic. So that the coefficients of Ti 
are e-periodic. Hence to get a total control of a coefficients of Ti, say ai, in a 
horizontal strip, we only need to control it in some vertical sub-strip on width £. 
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5.3.2. Analyticity and Fourier series. Moreover ai{ip) (and the other coefficients 
as well) is obviously analytic in by Theorem 3.1. This allows us to expand ai 
into an exponentially converging Fourier series. To control this series, we need to 
control both positive and negative Fourier coefficients, hence we need to have a 
control on the Wronskians in (5.14) et al. beyond the obvious zone Kq f] Ki (see 
Fig. 10). Actually we will get a uniform control in a strip of the form 

-A < Im < Im (^3 - (5, A > (5 > 0. 

Therefore we divide this strip into three different smaller strips called (I), (II) and 
(III) (see Fig. 11). Each of these strip will require a different type of computation. 

We will now start with describing the behavior of the main objects of our con- 
struction in each of the different strips. 

In the strip (II), we have a non-empty intersection for Kq and Ki. In the inter- 
section, we have 

(5.19) Ki = Kq, *i = -^^^ and = (j° . 

Consider the strip (I). Denote by 7 the common boundary of Kq and Ki in this 
strip (see Fig. 13). It is one of the Stokes lines beginning at = (^1. On 7, we 

get 

(5.20) «:i((^ + 0) = -«:o(¥^-0) 

(5.21) ^i(M,^ + 0) = i/^°(n,^-0) 

(5.22) (j4(<^ + 0) ='^°(<^-0) 

where 99 + (resp. 1^ — 0) denotes the limit taken from the right (resp. left). 
Formula (5.20) follows from (2.10), (5.21) from (2.6) and (5.22) from the definition 

of U)±. 

We also notice that the conformal properties of kq (see subsection 4.1) and (5.20) 
show that, for g(I) and Imi^ < 0, we have 

— ImKo(<^) = ImKi((^) < 0. 
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In the strip (III), the common boundary of Kq and Ki consists in the interval 
[(p2, 'Psl- This interval is the Stokes line joining (p2 to cp^. Along [(p2, 'Ps], we have 



5.3.3. A continuation lemma. To make the strategy outlined above work, we need 
to solve a technical problem that is that the asymptotics given by Theorem 3.1 
are only valid in admissible sub-domains of the canonical domains. So, if we want 
to continue the asymptotics across the boundary of a canonical domain to join, 
say, an admissible sub-domain of Kq to an admissible sub-domain of Ki, we will 
have to bridge a gap of width 2S. Therefore we will use the continuation lemma 
presented in the next subsection. We prove the 

Lemma 5.1. Let (^_, ip^ be two fixed points such that 

• Im = Im (^_|_, 

• there is no branching point of ip i^i'p) on the interval [<p-, <y5+]. 

Fix a continuous branch of k on [(p^,(p+]. Pick (po e ((/?_,(/?+). Let f{u,(p), 
f±{u,(p) be normal (i.e e-periodic in ip) solutions of (1-19) for ip e [(p^,(p+] and 
u e [—U,U] such that: 

1. f{u,(p) = e'^ ■^'fio'^'^^(^_^_(^u,(p) + o(l)) for (p e [(p^,(po] when e ^ 0, and the 
asymptotic is difjerentiable in u. 

2. f±{u,(p) = e^~' ^'^o'^'^''^{^j^(^u,Lp) + o{l)) for E[ip-,'P>j^] whene^O, and the 
asymptotic is differentiable in u. Here (^±) are chosen as in Theorem 3.1. 

Then, 

• if Im{K,{(p>)) > for all (p G [<P-, <p+], there exists C > such that, for e > 
small enough and (p G [(/?_, (p+], we have 



and the asymptotic is differentiable in u. 

Remark 5.1. This lemma is a precise version of a well known WKB heuristic 
that says that, if a function / admits on [(^_,(y9o] an asymptotic representation 
that is defined and exponentially increasing on [(^_,(^_|_], then the / admits the 
same representation on [(^o,(^+]. 

As throughout this paper, the estimates and asymptotics used for the solutions 
of equation (1.19) are differentiable in m, we will only write the relevant estimates 

df 

and asymptotics for the solution. We will not bother repeating them for — . 





• if Im{K,{(p)) < for all (p G [(^_, (^_|_], then, for (p G [(^_, (^+], we have 



f{u,^) = e^^"^.-''^{^^{u,^) + o{l)), 



du 
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As a last remark, let us say that in the rest of this section, to simplify the notations, 
we will denote by C every positive constant that is independent of (p and e. 

Proof. By assumption 2 of Lemma 5.1, (/-t) are linearly independent solutions 
of (1.19). So that we can write 



(5.27) 
where 



/ = a(^)/+ + %)/_ 



«(^) = :#4^ and %) "^^-'^^ 



^(/+,/-)' 



^(/+,/-) 

As in the proof of Theorem 3.1, we have 

\w{U, /_)| = |^>o|^=^o(l + = C ^ 0. 

Hence, on [(^_,(^o]) by assumptions 1. and 2., we have 

(5.28) a(^) = l + o(l). 

a is e-periodic in (p so that (5.28) holds on [(^_, (y^+j. 

Let us now estimate b. We start with the case when Im(K((^)) > for all 
if G [(^_,(^+]. Then, on [(^_,(^o]; we have 



h(/+,/)l<t^ 



g e JipQ ^ 



< 



So, for (/? G [<y3o — s,(fo], we get that \b{Lp)\ < C: hence, by e-periodicity, we have 
|^(<^)| < C for ip G [(^_, 99+]. On the other hand, on [(po, ip^], we have 

I/+I < Ce-^^^'o^'^W'^'^, |/_| < C'e^^^'o^'^W'^'^. 

So that, by (5.27), we get the announced result on [(^o,¥'+]- 

If lm{K{(pj) < for all (p G [<^_, (p+] then, for (p G [<^_, (po], we estimate 



|%)|<C 



Hence using this estimate for ip G [(^_ — e, (^_] and the e-periodicity, for ip G 

[^_,<^+], we get 

Hence, by (5.27), on [(^0;¥'+]; we have 
/ = a(^)/+ + %)/_ 
= geJcpo ^ I \l/+(n, (y^) + O ( e =■'''0 '"^j+o(l)l 

= e^^^'o'^'^^ (*+(«, ^) + 0(1)). 
This ends the proof of Lemma 5.1. □ 



48 



5.3.4. One more auxiliary canonical domain. To use Lemma 5.1, we introduce one 
more auxiliary canonical domain. We denote it by Kq. Let Kq be the continuation 
of Kq through the cut {ip2, (^3). Then Kq is defined to be the canonical domain for 
ko comprising the vertical interval {(p2, (fa). The part of Kq that will be needed to 
compute Ti is represented on Fig. 12. It will enable us to continue the asymptotics 
of a solution corresponding to the canonical domains Kq to the domain Ki along 
the horizontal lines limp = t] when Im (^2 < ^ < Im ■ 

We are now ready to start with the computations of the coefficients of Ti. 

5.3.5. The asymptotics of di{(p). We compute the asymptotics of w{f^, f]_) when 
£ ^ 0. 

Region (I). In view of (5.26) by Lemma 5.1, we 
can continue continue the asymptotics given 
by Theorem 3.1 for /° from [(^_,(/?o] to the 
whole interval [(^_.(^+] (see Fig. 13). The 
consistent basis needed to apply Lemma 5.1 is 
obtained by noticing that, as /«o(^) = —i^oi'f), 
the domain Kq symmetric of Kq with respect 
to the real axis also is a canonical domain. 
This domain comprises a neighborhood of 7, 
the common boundary of Kq and Ki in region (I) (see Fig. 13). 
By what has been said above, on [ipi,, (p^], we have 

By (5.17) and as Ki = —Kq in the strip (I), we know that 

/_j_ = e=-''^ ^ =■'''1 ^■'^b ^[^^ + o{l)). 

The curve 7 joining ipi to (ph is a Stokes line for kq (and ni), hence 

Moreover we notice that for G [(^i,7r], Ko{ip) e R. 

Now as if has to be in an admissible sub-domain of Kq, we know that Re (ip — iph) > 
S; hence, we get that 

(5.29) \w{f_,fl)\ < C\e-^^o-od^\eCS/^ = c\e-^ ''°'^^\e^^^' . 

This was done on the segment [(pb,^+]- By 
e-periodicity, it will hold everywhere on the 
line {Im((^) =lm((ph)}. The estimate obtained 
above is valid as long as the point (p^ stays 
in the admissible sub-domain of Kq and if 
the interval [(^_ , ip+] stays off some vicinity 
of the branching points. Hence, we get that, 
for —A <luiip< —S (A, c) > 0), we have 

(5.30) \di{(p) \ < C|e-^^o'''^°'^^|e^^/^ 
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Region (II). In this case, by (5.19), for ip & KqCi Ki, we have 
Comparing this with the formula 



i/o"«o^^e^o"'"-'^^(u;(/^, /i) + o(l)). 



one easily obtains 

Hence, as long as we stay away from the turning points (pi and (p2, in region (II), 
we get that 

(5.31) di{^) = e-i-^o''°'^^e-fo''-'^^{l + o{l)), S <lm(p <lm^2 - ^ (S > 0) 
Here the integrals are taken along a curve in KqC] Ki. 

Region (III). In this region, we will satisfy 
ourselves with an estimate on w{f°, fl_). As 
Im/^o > in this sector, by Lemma 5.1, we 
can continue the asymptotics given by The- 
orem 3.1 for /° from [(p_,(^o] to the whole (II) 
interval [(/?_, (^+] (see figure 14). The consis- 
tent basis needed to apply Lemma 5.1 is ob- 
tained by using Theorem 3.1 for the canonical 
domain Kq. This domain comprises a neigh- 
borhood of the common boundary of Kq and Ki that is the segment ((^2, ^3) (see 
Fig. 11, 12 and 14). On the segment [(pi,,(p+], we have 

Hence, for ip G (^+], we get 

\f-\<C\e- 

We also have 




\fl\<c 



To get this estimate, we have used (5.23). The last two estimates imply 



e ^ 



\dM\<C 

Now we note that 

• for (y^i < (y3 < TT, /Co e M. 

• for ip e [(p2, fs], Ko = + 'it where t > 0. 
Assuming that \ip — iph\ is of order 5, we obtain 

(5.32) Mi(<^)| < C|e-^^o'''^° 

Since di is e-periodic, this estimate is valid 
in the strip {lm(p2 + S < limp < Im (p^ — S} 
{S>0). 
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Conclusion: as di{(p) is e-periodic in (p and analytic in a strip {—A < limp < 
Im(^3 — 6}, we can expand di in a Fourier series with exponentially decreasing 
coefficients 

where, for any po G {—A < Imp <lmps — S}, 

Sn = - d,{p)e-^'^^^dp. 

^ Jipo 

Using estimates (5.30) and (5.32) in the regions (I) and (III), we get that, for 
n ^ 0, we have 

\Sn\< C\e-^ -^0 ''°'^V\f,-2n\n\Yo/e^C6/e _ 

where Yq = min(A, Im p^ — 5, Y) and 5 is chosen such that C5 < 27r(F — Yq). 
Using estimate (5.31) in region (II), we obtain 

Then, uniformly for p in the strip {|Im2;| < Y}, we have 

di{p) = e-?/o"'=od^+/o"^°'iv(i + o(l)) 

5.3.6. The asymptotics of bi{p). We will only estimate w{f1, /|) when e ^ 0. 

Region (I). Let 7 be the Stokes line beginning at p) = pi and bordering Kq and 
Ki. As Imno > 0, Lemma 5.1 will only give us an estimation on when we cross 
7 along [(^_, p+]. This estimation is 

As on 7 and on [p_, p+], kq = —Ki, we can estimate by 

I/+I < C|e^^''"°^^"^^''"°^^"^C"°'^^||e"^^^'o"<"^^|. 
These two estimates imply 

|^'^(/+,/+)| < C\e^^o''^odv\\e^i:^\imno\d^^_ 

for p on [ph, p+] and inside the admissible domain for Ki. Now the only restriction 

we have on po is that it has to be in the admissible domain for Kq] hence we know 
that we can choose |</7o — ^| ^ ^ where S is the distance that defines the admissible 
domains (see (3.5)). It may be chosen as small as needed. So, for — A < Im p < —5 
(A, ^ > 0), we have 

(5.33) \bi{p)\<Ce'^^ei^o''od^\. 

Region (II). As we only estimate Bi{p), we don't need to compute it in the region. 

Region (III). As Im(Ko) > 0, Lemma 5.1 will only give us an estimation on in 
Ki. The point p)b being as on Fig. 14 the intersection of [p2, pa] and [p^, p+], we 
get 
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for being in an (5-vicinity of lpq. Similarly, we have 
under the same conditions on ip. 

We know that along [p>2, ¥^3] and on [tpo, we have kq = 271 — ki] hence 

\w{fl, fl) I < C|ee "o^^^l |e'^ S^^ ^^+2-^ C '^^'i |gf I^^ \imno\d^^ 

as (p — (fo > S . We now use the facts that kq is real positive along [tt, ^2] and that 
Im((^h — P2) > 0. Hence, in (III), we have 

|6i(^)| < Ce^ilee^o^^o^^le-"/" 

where 77 is a fixed positive constant (independent of both S and s) as Im (JJ^ koc/<^+ 
J^^ dp) > Hence, for S small enough, 

Conclusion: using the same method as before, in the strip {|Im(^| < Yq}, we get 
that 

5.3.7. The asymptotics of ai{p). We compute the asymptotics of w{f^, /I) when 
£ ^ 0. The computations are done in the same way as for bi and di. 

Region (I). We get that, 

\aM \ < Ce^^e^-^o'^od^l, 
for -A <lmp < -S (A, ^ > 0). 
Region (H). We get that 

ai(^) = ei^o''°'^^e^o<'^^{l + o{l)), 
for < S < Im p < Im P2 — S . 
Region (HI). We get that, 

for Im + < Im < Im 993 — 5 (^ > 0). 

Conclusion: estimating the Fourier coefficients of ai, as above, we get that, uni- 
formly for p in the strip {|Im2;| < Y}, 

a,{p) = e^/o"«'"^^+/o"-V^(l + 0(1)). 

5.3.8. The asymptotics of Ci{p). We compute the asymptotics of w{f'^,f^) when 
e ^ 0. 

Region (I). As Im(«:o) > 0, we can use Lemma 5.1 to continue the asymptotic 
expansion for /° given by Theorem 3.1 in Kq across OKq along a horizontal line 
(see Fig. 13). This will allow us to get a precise asymptotic of Ci. We use formulae 
(5.20) - (5.22). The asymptotic expansion for /° on [(^6, (^+] is 

f_{u,p) = e-^^o^'^od^ . (^yA;^e/o"'-°'i^(^o +0(1))) . 
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Here we integrate along paths going from to (p staying below the point (pi. We 
use (5.20) and compute 



(5.34) 



pip r-ipi r-lf 

I Kodip = / Kodip — / Kidip 

Jo Jo Jipi 



/ Kodip — / Kidip — / Kidip. 

JO J fi J IT 



In the first two integrals , we integrate along the real axis; in the last integral, the 
contour is located in Ki. Similarly 

pip pifi pifi 

/ co'^dip = / uP_dip + / uj\d(f, 

Jo Jo Jifi 

where the first integral is taken along the real line; the integration contour in the 
second integral lies in Ki. Finally, we get 

pifi pifii pipi pip 

(5.35) / uj'^dip = / uP_d(f — I uj\dip -\- / co^dip, 

Jo Jo Jn Jn 

We also note that, on the difl'erent sides of 7 in region (I) (see Fig. 13), we have 

(5-36) V^i^+o = '^y/^iip-o 

Combining (5.34) - (5.36), we get 

fO/ \ ■ fif^ Kodip+^ r Kodifi+fif^ 0J° -1"^^ L0° i f«'Kidv/,T,W N , /tw 

j_[u,(p) = —ie ^ ■'° ''^^BJipi ^ ^^jo - ^^["^^[u,!^) + o[l)). 

On the other hand, we know that 

Hence we easily obtain 

that is 

(5.37) ci(v?) = -i . e--^ ^0 "^^^e^i ^^i''°'^^e^°'' + o{l)) 

where the integrals are taken along the real line. 

This holds in the strip {-A <lm(p< -S} {A > S > 0). 

Region (H). As for bi, we will not need any information on Ci in this region. 

Region (HI). As Im(/co) > 0, we will again be able to use Lemma 5.1 to compute 
a precise asymptotic for ci. So that, for ip G [(p-, (^+], we know that 

(5.38) f_ = e-^ ^o^'^orf^ (v^e^o"'^°'i^^o ^ ^^^^^ 

where [(^_,(^_|_] is defined in Fig. 14. In (5.38), the integration contour first goes 
from to (p_ in Kq, then to ip along [(^_,(^+]. Using (5.24) - (5.25), we rewrite 
the integrals in (5.38) to get 

pip pir p(p2 pip 

/ Kodip = / Kodip + 2 / acqc/i^ + 27r((^ — (^2) — / nidip, 

Jo Jo Jtt Jit 

where, in the first integral, we integrate along a path in KqCiKi, and in the next 
two, along a path in Ki. 
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The same way we get 

^0 Jo Jtt Jtt 

We also note that, for 99 e [(^2, ^3], we have 
(5.39) 

Combining these formulae, we obtain 



Using this in conjunction with 
we get 

(5.40) ci((^) = ie'Uo ''°'^'P-mpif^o-7^)d'P^-2niiip-n)^jr ujUf-Jp ujldv(^l _^ ^(^i-^y 

This asymptotic is valid for Im <^2 + < <y5 < Im </?3 — 5 {S > 0) . 

Conclusion: now to get global information on Ci, we use the same method as 
above; let 



Cl 



Then, by the asymptotic (5.37), for n > 0, we have 
The same formula gives 

= ^^-Uo -odv+2i Kod^g/o"! u.U'p-j:^ ^id^^i + 0(1)). 

We use the asymptotic given by region (III) to estimate (7^) for n < 0. For some 
(po in (HI); we compute 



1 f* 

7-1 = - / ci(^)e^^'^"= dip 

= i^-Uo >^odv-H SPi>^o-^)dv^Sr u;Uv-SP c^ld^f^i + 0(1)) 



Now we recall that the segment [(^2,<^3] is a Stokes line for kq and ko(<^2) = tt] 
hence Im J^^{ko — 7r)d(p = 0. For n < —1, we use (5.40) to estimate 



1 7" 



£ / 



ci{ip)e^'''^ d<p 



< C 



/ ' 

JtpQ 
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Putting all this together, we get that, for ip in the strip {|Im2;| < Yq}, we have 



where o(l) is uniform in the strip {|Im2;| < Yq\. This ends the proof of Proposi- 
tion 5.1. 

■5.4. The asymptotics of T2 . As we have already written, we will not do the 
detailed computations for T2. We will only describe the geometry that is slightly 
different. Let us write 



C2W) d2{^))' 



then 
(5.41) 

(5.42) 

(5.43) 



a2{^) = -w{fl /!), b2{^) = --w{fl fl), 

W2 W2 

C2{^) = -W{f_,f_), d2{^) = --Wiflfl), 
W2 1^2 



T2 shares many of the properties of Ti (though not in the same domain); in 
particular: 

• its coefficients are e-periodic, 

• they are analytic in some strip. 

Now due to the geometry of Ki and K2 (see Fig. 15), we will only be able to 
control these coefficients in some strip of the form 

-lm(ps + S <lm(p < A, A>S>0. 




Again we divide the strip {—Imip^ + S <lianp < A} into three smaller strips 
called respectively (I), (II) and (III) (see Fig. 15). 

Let us briefly describe the main objects of our construction in each of these strips. 

In the strip (II), we have a non-empty intersection for Ki and K2. In the inter- 
section, we have 

1^1 = 1^2, '01 = '0± and Lo]. = uj\. 

Consider strip (I) . Denote by 7 the common boundary of K2 and Ki in this strip 
(see Fig. 16). It is one of the Stokes lines beginning at = 27r — 991. On 7, we get 

Ki{ip - 0) = -K2('y5 + 0) 

ii{u,ip-^ = ni^,^ + {)) 
a;i(^-0) =0)2(^ + 0). 

These formulae immediately follow from equa- 
tions (5.20) - (5.22) as, for (p G Kq, we have 

(5.44) K2{2'K - if) = Kq{^). 

We also notice that the conformal properties 
of Kq (see subsection 4.1) and (5.20) show 
that, for g(I) and Im(^ < 0, we have 

— ImK2(<^) = ImKi((^) < 0. 

In the strip (III), the common boundary of K2 
and Ki consists in the interval [27r — (^2, 27r — 
Lp^]. This interval is the Stokes line joining 
27r — Lp2 to 27r — (^3. Along [27r — (^2, 27r — (^3], 
we have 

Ki(^-O) = 27r- ^2(^ + 0) 

^i(?i,^-0) = ^J(m,^ + 0) 

a;l(^-0) =a;^(^ + 0). 

As in region (I), for — Im (^3 < Im < — Im ip2, 
we have 

— ImK2((^) = ImKi((^) < 0. 

To use Lemma 5.1 to compute T2, we will use another auxiliary canonical domain, 
namely the symmetric of Kq with respect to the point tt. It comprises a neigh- 
borhood of (27r — (^2, 27r — (^3). It is canonical with respect to K2 the branch of the 
quasi-momentum obtained from k,2 by analytic continuation through the interval 
[27r — ^p2^ 27r — (^3] by (5.44). How to use Lemma 5.1 is partly described in Fig. 16 
and 17. We also point out that, in this case, (5.39) becomes 




We will not detail the computations of the coefficients of the matrix T2 any further. 
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Remark 5.2. Let us just notice that one could have used the symmetry between 
the domain Kq U Ki and KiU K2, and between the associated consistent basis to 
compute T2 from Ti. 

5.5. The computation of M. Let K := Kq U Ki U K2. This domain is rep- 
resented on Fig. 18. Its natural cuts are the branches of Stokes lines shown on 
Fig. 18 by bold lines. On this domain, we define a unique branch of the complex 
quasi-momentum k by 

This branch is analytic on K. In the same way, on K, we define a unique analytic 
function co by 

co± := co^ = ujI_. 




These branches have the following symmetry properties: 

(0.45) K{27r — (fi) = k{(p) and iO±{27r — ip) = —uj±{(p). 

The first equality is just (4.9); the second is obtained using (3.6) and the fact 
that, as cos(27r — ip) = cos{(p), we have p±{x, 2% — ip) = p±{x, (p). 

Rewriting the asymptotics of the matrices Ti and T2 in terms of k and u^. and 
using the symmetry property (5.45), one computes the product of these matrices. 
It is an elementary but long and tedious calculation that we will not reproduce 
here. This ends the proof of Theorem 5.1. □ 

We end this section with proving that 02 is real. Therefore we compute 
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I-2-K — IP2 

as u)J^{^p) = ijj^ilp), ^2 = 27r — 1^2 and / (a;+ — u)J) = 0. 

6. Spectral application 

6.1. Reduction to a difference equation. Here, we discuss the passage from 
the family of equations (1.1) to the difference equation (1.9) and justify the relation 
(1.11). 

Consider the vector \& composed of the consistent basis solutions of (1.1) as in 
(1.7). The definition of the monodromy matrix implies that 

* {x+2Tm/£, (j) + 2Tm/e) = 

(6.1) = M ((/) + 27r(n - l)/e) M (</) + 27r(n - 2)/c) • • • 

•••M(0 + 2^/c)M((f))*(a-,</)). 

Assume that x is a fundamental solution of equation (1.9) defined by the mo- 
nodromy matrix M corresponding to By definition, x(</) + h) X~^(0) = M[(j)). 
Since M{(l) + 27rl/£) = M{4) + lh) (for I e N), together with (6.1), this implies the 
relation 

* (x + 27rn/£, (f) + 27rn/£) = x{(J) + nh) x'^{4>) * {x, cj)). 
Since is 1-periodic in 0, it immediately leads to (1.11). 

6.2. Bloch solutions of difference equations. 

6.2.1. To prove Theorem 1.3 we have to construct some solutions of equation 
(1.9). So, we recall some basic ideas (see [8]) related to difference equations of the 
form 

(6.2) 4){^ + h) = M{(l))4){^), (/)GM, 

where /i is a positive number and M {(/)) e SL (2, C) is a given matrix function 
1-periodic in </). 

U ip is a fundamental solution (i.e. detip = 1), then any other solution can be 
uniquely represented in the form 

x(</.) = ^((/.)P (</.), <^eM. 

where P is an /i-periodic matrix function, P {(f) + h) = P {(j))- In particular, one 
can write 

(6.3) ^){cj)^l) = ^j{(i))Ml {(fy. 

This relation defines the monodromy matrix Mi corresponding to the fundamental 
solution ip. Clearly, 

(6.4) Mi((/) + /i) = Ml ((/.), det Ml ((/)) = 1, V</) e E. 

If the monodromy matrix corresponding to a fundamental solution is diagonal, 

^^^'^> - \ (i-\cj>) ) ' 
we call tp, a Bloch solution, and the coefficient /5 is its Floquet multiplier. 
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We pick the determination of the function In (5 {(f)) so that all its possible discon- 
tinuity jumps be smaller than 27r. If the limit 

1 " 

(6.5) U) = lim - V In ^ (</) + I) 

1=0 

exists and its convergence is uniform in (/), we call this limit the complex Lyapunov 
exponent of the Bloch solution tp. 

We call a Bloch solution monotonous if Re 9 {(f)) > Const > for all (f). 

We finish this section by formulating a simple condition for the existence of the 
monotonous Bloch solutions [8]. 

Set 

(6-6) P('^) = Tr%^' v{(f>) = Mn{cf))+p{c^)M,,{(f)-h). 

Mi2(0 - h) 

Let 

inf|p(</')|, p+ = sup v^ = ml\v{(f))\. 



For / (0) 7^ being a periodic continuous function on R, we denote by ind / the 
integer equal to the increment of over one period. Clearly, ind/ is independent 
of (f>. One has 



Proposition 6.1. Let M be continuous in (f), and let M12 ^ 0. If 

2 

p 



> 0, ^ ("2") ' ind f = ind p = 0, 

then 

1. equation (6.2) has a fundamental Bloch solution continuous in (f); 

2. the Lyapunov exponent satisfi,es the estimates 

ln(«_/2 + V(^-/2)2 - p+) <hRe{e {(f))) < \n{v+ + i;_/2 - V(^_/2)2 - 

6.3. Bloch vectors for the family of equations(l.l). In terms of Bloch so- 
lutions of (1.9), one can construct a very natural set of solutions of the family of 
equation (1.1). Let be a vector composed of consistent basis solutions of (1.1), 
and let M be the corresponding monodromy matrix. Suppose that x is a Bloch 
solution of equation (1.9). Set 

F{x, (f>)=x-\<f>)'^{x,<f>). 

The components Fi and F2 of the vector F satisfy (1.1). First, we note that 
w (Fi,F2) = detx~^{(t)) w (*i, ^^2), and thus, 

(6.7) w{FuF2) = l, V</)eM. 
Written in terms of F, formula (1.11) takes the form 

(6.8) F{x + 27rn/£,(f)) = F{x,(f)-nh), Vx, (/) e M. 
Finally, one can easily see that 

(6.9) F{x,(f) + 1)= ^"q^^^ ^ F {x, </.), Vc/. G M, 
where /5 is the Floquet multiplier of the Bloch solution x- 
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We shall call any vector of solutions of (1.1) possessing the properties (6.7)-(6.9) 
a Block vector of this family of equation, and the coefficient 1//3 from (6.9) its 
Floquet multiplier. 

6.4. Bloch vectors and the spectra of (l-l). As it was done above for Bloch 
solutions of equation (6.2), we define the complex Lyapunov exponent of a Bloch 
vector of (1.1) in terms of its Floquet multiplier. If the Lyapimov exponent verifies 
the estimate of the form Re^ (x) > Const > for all x, we call the Bloch vector 
monotonous. 

Clearly, if there is a monotonous Bloch solution of (1.9), then there exists a mo- 
notonous Bloch vector of (1.1). We shall use the following elementary observation. 

Lemma 6.1. Assume that, for a given E, there exists a monotonous Bloch vector 
for (1.1); n,nd that its components are locally bounded in x and (j). Then, for any 
(/) e M, E does not belong to the spectrum of (1.1) in L^(M). 

Proof. For any given (j), Fi^2 are solutions of (1.1). By (6.7), these solutions are 
linearly independent. Estimate the first of them. In view of (6.8) - (6.9), 

N-l 



\ 1=1 / 

Since the limit (6.5) exists and is uniform in x, and the Bloch vector is monotonous, 
then 



for sufficiently large n with some positive S. In the same way, one checks that 



for sufficiently large — n. The linear independence of the solutions Fi 2 and their 
estimates imply that the energy E is in the resolvent set of the operator defined 
by the left hand side of (1.1) (see [10]). 

6.5. Exponential localization of the spectrum. Now, we can immediately 
prove Theorem 1.3. Consider equation (6.2) with the matrix M described in the 
Theorem 1.2. In this case, 

p = 1 + O (1), V (z) = \ (C0S((/.1 + O (1)) - ti (1 + O (1)) COS Z ) (1 + (1)) 

with z chosen as in Theorem 1.2. All the error terms are locally analytic in E. 
The assumptions of Proposition 6.1 are satisfied if 



This condition guarantees the existence of monotonous Bloch solutions of (1.9). 
Now, combining Proposition 6.1 and Lemma 6.1, we see that, for E being in the 
interval described by (1.26), the spectrum of (1.1) can be situated only on the 
subintervals where 



Fi(3; + 2n7r/e,(/))| = ]^ |/5 (</) + /)||Fi(a;, - n/i + 7V)|. 



Let < X < 27r/e. Choose N so that 0<</)-n/i + Ar<l. Then 




Fi(a; + 2n7r/£,0)| < Ce 



F2(x + 2n7r/e,</))| < Ce 



cos((/.i + o(l))| > (t + ti) (1 + 0(1)). 



(6.10) 



cos(</.i + o(l))| < (t + ti) (1 + 0(1)). 
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Remind that t, ti are analytic in E' in a vicinity of any point E satisfying 
(1.26). We show that 

(6.11) <f>'i{E)^0, 

for any E satisfying (1.26). This function is defined by (5.1). It can be represented 
in the form 




Ki dep. 



Here Hi is the branch of the complex momentum corresponding to the canonical 
domain Ki, see section 5. Recall that it is related to a branch of the Bloch 
quasi-momentum of (1.14) by the formula 

Ki {(fi) = k {E — cos (fi) . 

So that, for E satisfying (1.26), 

(6.12) (l)'AE) = - / k' (E - cos(p)d(p. 

Omitting the elementary details, we only note that, since A; as a function of ip 
has square root branching points at ip = cpi and cp = 27r — ipi, this integral (6.12) 
converges. Moreover, clearly, up to its sign, it does not depend on the choice of 
the branch of k. For any ip G {(pi, 27r — ipi), 

El < E — cos (p. 

Moreover, under the condition (1.26), 

E - cosip < E2 - S, 6 > 0, ip eR 

As k'{E — cosip) is real and these two inequalities imply that k'{E — cos(p) ^ 
on the interval {ipi, 27r — ipi), we get the desired result. 

The proof of Theorem 1.3 follows from (6.10) - (6.11). We now establish the 
properties of the subintervals described by (6.10). Since t and ti are exponentially 
small and (p[ / 0, the length of these intervals are exponentially small as well. 
They are situated in o (£:)-vicinities of the points E^'''> defined by 

cos(/.i(E(^)) = 0. 

Since (f)[ 7^ 0, then the distances between these points are of order e, and near 
each of these points, there is only one of the subintervals containing the spectrum. 
Recalling the definitions of t and ti (cf (5.1)), one sees that, up to the factor 
1 + 0(1), these functions are constant on the exponentially small intervals. This 
implies that the length of the subinterval situated near to E^''^ is equal to {t{E^''^) + 
ti(£'(^^)) (1+0(1)). This completes the proof. □ 
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